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HARDY SPACE THEORY ON SPACES OF HOMOGENEOUS TYPE 
VIA ORTHONORMAL WAVELET BASES 


YONGSHENG HAN, JI LI AND LESLEY A. WARD 

Abstract. In this paper, using the remarkable orthonormal wavelet basis constructed re¬ 
cently by Auscher and Hytonen, we establish the theory of product Hardy spaces on spaces 
X = X\ x X 2 x • • • x X n , where each factor X t is a space of homogeneous type in the sense 
of Coifman and Weiss. The main tool we develop is the Littlewood-Paley theory on X, 
which in turn is a consequence of a corresponding theory on each factor space. We define 
the square function for this theory in terms of the wavelet coefficients. The Hardy space 
theory developed in this paper includes product H p , the dual CMO p of H p with the special 
case BMO = CMO 1 , and the predual VMO of H 1 . We also use the wavelet expansion to es¬ 
tablish the Calderon-Zygmund decomposition for product H p , and deduce an interpolation 
theorem. We make no additional assumptions on the quasi-metric or the doubling measure 
for each factor space, and thus we extend to the full generality of product spaces of homo¬ 
geneous type the aspects of both one-parameter and multiparameter theory involving the 
Littlewood-Paley theory and function spaces. Moreover, our methods would be expected to 
be a powerful tool for developing wavelet analysis on spaces of homogeneous type. 
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1. Introduction 

We work on wavelet analysis in the setting of product spaces of homogeneous type in the 
sense of Coifman and Weiss |CWlj . where each factor is of the form (A", d, p) with d a quasi¬ 
metric and p a doubling measure. We make no additional assumptions on d or p. After 
recalling the systems of dyadic cubes of Hytonen and Kairema [ HK ] and the orthonormal 
wavelet basis of Auscher and Hytonen [AH], we define an appropriate class of test functions 
and the induced class of distributions on spaces of homogeneous type. We prove that the 
Auscher-Hytonen wavelets are test functions, and that the Auscher-Hytonen reproducing 
formula for L p also holds for our test functions and distributions. We show that the kernels of 
certain wavelet operators Z\, defined in terms of these wavelets satisfy decay and smoothness 
conditions similar to those of our test functions. These facts play a crucial role in our 
development of the Littlewood-Paley theory and function spaces, later in our paper. 

We define the discrete Littlewood-Paley square function via the Auscher-Hytonen wavelet 
coefficients. In order to establish its L p -boundedness, we also introduce a different, continuous 
Littlewood-Paley square function defined in terms of the wavelet operators D^. We prove 
that the discrete and continuous square functions have equivalent norms, by first establishing 
some inequalities of Plancherel-Polya type. We develop this Littlewood-Paley theory first in 
the one-parameter setting, and then for product spaces. 

For p in a range that depends on the upper dimensions of the spaces X\ and X 2 and strictly 
includes the range 1 < p < 00 , we define the product Hardy space H p (X\ x X 2 ) as the class 
of distributions whose discrete Littlewood-Paley square functions are in L p (Xi x X 2 ). (Here 
we write only two factors, for simplicity, but our results extend to n factors.) For p in this 
range with p < 1, we define the Carleson measure space CMO p (Xi x X 2 ) via the Auscher- 
Hytonen wavelet coefficients, as a subset of our space of distributions, and prove the duality 
(H p (X\ x X 2 ))' = CMO p (Xi x X 2 ) by means of sequence spaces that form discrete analogues 
of these spaces. This duality result includes the special case (H 1 (Xi x X 2 ))' = BMO(Xi x X 2 ). 
We define the space VMO(Xi x X 2 ) of functions of vanishing mean oscillation, also in terms 
of the Auscher-Hytonen wavelet coefficients, and prove the duality (VMO(Ad x X 2 ))' = 
H 1 (X 1 x X 2 ) by adapting an argument of Lacey-Terwillcger-Wick [LTW] . Using the wavelet 
expansion, we establish the Calderon-Zygmund decomposition for functions in our Hardy 
spaces H p (Xi x X 2 ), again for a suitable range of p that strictly includes 1 < p < 00 . As 
a consequence, we deduce an interpolation theorem for linear operators from these product 
Hardy spaces to Lebesgue spaces on X\ x X 2 . 

We now set our work in context. As Meyer remarked in his preface to [DEI, “One is 
amazed by the dramatic changes that occurred in analysis during the twentieth century. In the 
1930s complex methods and Fourier series played a seminal role. After many improvements, 
mostly achieved by the Calderon-Zygmund school, the action takes place today on spaces 
of homogeneous type. No group structure is available, the Fourier transform is missing, 
but a version of harmonic analysis is still present. Indeed the geometry is conducting the 
analysis. ” Spaces of homogeneous type were introduced by Coifman and Weiss in the early 
1970s, in [CW1] . We say that (A", d, p) is a space of homogeneous type in the sense of Coifman 
and Weiss if d is a quasi-metric on X and p is a nonzero measure satisfying the doubling 
condition. A quasi-metric d on a set A" is a function d : X x X —* [0, 00 ) satisfying (i) 
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d(x,y ) = d(y,x ) > 0 for all x, y G X; (ii) d(x,y ) = 0 if and only if a; = y\ and (iii) the 
quasi-triangle inequality, there is a constant A 0 G [1, oo) such that for all x, y, z G A", 

(1.1) d(x,r/) < A 0 [d(x,z) +d(z,y)]. 

We define the quasi-metric ball by B(x,r ) := {y G X : d(x,y) < r } for x G X and 
r > 0. Note that the quasi-metric, in contrast to a metric, may not be Holder regular and 
quasi-metric balls may not be open. We say that a nonzero measure p satisfies the doubling 
condition if there is a constant C M such that for all x G X and r > 0, 

(1.2) p(B(x,2r)) < C fl p(B(x,r)) < oo. 

We point out that the doubling condition (11.21) implies that there exist positive constants to 
(the upper dimension of p) and C such that for all x G X, A > 1 and r > 0, 

(1.3) p(B(x, Xr)) < CX ul p(B(x, r)). 

Spaces of homogeneous type include many special spaces in analysis and have many appli¬ 
cations in the theory of singular integrals and function spaces; see [CW21 INSll 1NS2] for more 
detail. For instance, Coifman and Weiss introduced the atomic Hardy space on (A, d, p) and 
proved that if T is a Calderon-Zygmund singular integral operator and is bounded on L 2 (X), 
then T is bounded from H P (X) to L p (X) for some p < 1. 

However, for some applications, additional assumptions were imposed on these general 
spaces of homogeneous type, because as noted above the original quasi-metric d may have 
no regularity and quasi-metric balls, even Borel sets, may not be open. For example, to 
establish the maximal function characterization of the Hardy space introduced by Coifman 
and Weiss, Macias and Segovia proved in [MS} that one can replace the quasi-metric d by 
another quasi-metric d! on X such that the topologies induced on X by d and d! coincide, 
and d! has the following regularity property: 

(1.4) \d'(x,y) - d'(x',y)\ < C 0 d'(x, x') e [d'(x, y) + d'(x', y)} 1 ^ 

for some constant C 0 , some regularity exponent 6 G (0,1), and for all x, x', y G X. Moreover, 
if quasi-metric balls are defined by this new quasi-metric d', that is, B’(x,r ) := {y G X : 
d'(x,y) < r} for r > 0, then the measure p satisfies the following property: 

(1.5) p(B'(x,r)) ~ r. 

Note that property (1 1.5ft is much stronger than the doubling condition. Macias and Segovia 
established the maximal function characterization for Hardy spaces H p (X) with (1 + $) _1 < 
p < 1, on spaces of homogeneous type (A", d', p) that satisfy the regularity condition (j 1,4ft on 
the metric d! and property (jl .5[) on the measure p\ see [MS) . 

A fundamental result for these spaces (A ,d’,p) is the T{b) theorem of David-Journe- 
Sennnes [DJSj . The crucial tool in the proof of the T{b ) theorem is the existence of a 
suitable approximation to the identity. The construction of such an approximation to the 
identity is due to Coifman. More precisely, take a smooth function h defined on [0, oo), equal 
to 1 on [1,2], and equal to 0 on [0,1/2] and on [4, oo). Let X/ be the operator with kernel 
2 k h(2 k d\x,y)). The property (j 1.5 [) of the measure p implies that C~ l < T/(l) < C for some 
C with 0 < C < oo. Let M k and W k be the operators of multiplication by 1/T/(1) and 
{T fc [l/T fc (l)]} _1 , respectively, and let S k M k T k W k T k M k . It is clear that the regularity 
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property (II,4j) on the metric d! and property (1 1.5ft on the measure /i imply that the kernel 
S k (x,y ) of Sk satisfies the following conditions: for some constants C > 0 and £ > 0, 

(i) S k (x,y) = 0 for d'(x,y) > C2~ k , and H^Hoo < C2 k , 

(ii) | S k (x,y) - S k (x',y)\ < C2 k(1+£) d'(x,x'Y, 

(ih) \Sk(x,y) - S k (x,y')\ < C2 k(1+e) d'(y, y') £ , and 

(iv) / S k (x,y)dfj,(y) = 1= / S k (x,y) dy{x). 

J x J x 

Let D k := S k+1 - S k . In (EE], the Littlewood-Paley theory for L P (X), 1 < p < oo, was 
established; namely, if y(X) = oo and y(B(x,r)) > 0 for all i£l and r > 0, then for each 
p with 1 < p < oo there exists a positive constant C p such that 

cy'n/iip < l|{Ei°‘(/)i 2 hllp^Gii/iip- 

k 

The above estimates were the key tool for proving the T(l) theorem on (A", d', /i); see [DJSl | for 
more detail. Later, the Calderon reproducing formula, test function spaces and distributions, 
the Littlewood-Paley theory, and function spaces on (A ,d',p) were developed in |Hlj . (HS ] 
and (H2). However, in those works wavelet bases were replaced by frames, which in many 
applications offer the same service; see ra for more details. 

In (NSlj . Nagel and Stein developed the product L p (1 < p < oo ) theory in the setting of 
the Carnot-Caratheodory spaces formed by vector fields satisfying Hormander’s finite rank 
condition. The Carnot-Caratheodory spaces studied in [NSlj are spaces of homogeneous 
type with a smooth quasi-metric d and a measure y satisfying the conditions y(B(x, sr)) ~ 
s m+2 y(B(x, r)) for s > 1 and y(B(x, sr)) ~ s 4 y(B(x,r)) for s < 1. These conditions on 
the measure are weaker than property (11.5ft but are still stronger than the original doubling 
condition (11.2ft . In [ HMYj . motivated by the work of Nagel and Stein, Hardy spaces were 
developed on spaces of homogeneous type with a regular quasi-metric and a measure satisfying 
the above conditions. Recently, in HLL2I , it was observed that Coifman’s construction of an 
approximation to the identity still works on spaces of homogeneous type (A, d, p) where the 
quasi-metric d satisfies the Holder regularity property (jl . 4(1 but the measure // only needs to 
be doubling. Specifically, the kernel S k (x,y) of the approximation to the identity S k satisfies 
the following conditions: there exist constants C > 0 and 6 > 0 such that for all k G Z and 
all x, x ', y, y f G A, 


(i) S k (x,y ) = 0 for d(x,y) > C2 fe , and |S' fc (x,i/)| < C 


V 2 - k {x) + V 2 - k {yY 


v 


j\0 


(ii) \S k (x,y) - S k (x',y)\ < C2 ’ d{x,x') .. , 

V 2 -k{x) + V 2 -k(y) 

(iii) property (ii) also holds with x and y interchanged, 

(iv) |[Sfc(a;, 2 /) - S k (x,y')\ - [S k (x',y) - S k (x',y')]\ 

<C2 2 “d(i,x')''d(a,a')‘ i 1 


v 2 -,(x) + v 2 My)’ 

Sk(x,y) dy(y) = 1 = / S k (x,y)dy(x), 


and 


IX 


where V r (x) := y(B(x,r)). 

































HARDY SPACE THEORY ON SPACES OF HOMOGENEOUS TYPE VIA WAVELET BASES 


5 


For spaces of homogeneous type (A", d, /i) with some additional assumptions, the one- 
parameter and product Hardy spaces were developed in [HMYj and 1111.1/2 . respectively. 

A natural question arises: can one develop the theory of the spaces H p and BMO on spaces 
of homogeneous type in the sense of Coifman and Weiss, with only the original quasi-metric d 
and a doubling measure /i? 

Recently, Auscher and Hytonen constructed an orthonormal wavelet basis with Holder reg¬ 
ularity and exponential decay for spaces of homogeneous type [SHI- This result is remarkable 
since there are no additional assumptions other than those defining spaces of homogeneous 
type in the sense of Coifman and Weiss. 

Auscher and Hytonen’s orthonormal wavelet bases open the door for developing wavelet 
analysis on spaces of homogeneous type in the sense of Coifman and Weiss. Motivated 
by Auscher and Hytonen’s work, the purpose of the current paper is to answer the above 
question. Specifically, we will employ a unified approach to establish a product Hardy space 
theory on X = X\ x X 2 x • • • x X n , where each factor is a space of homogeneous type in the 
sense of Coifman and Weiss. It was well known that any analysis of the product Hardy space 
on X — Xi x X 2 x ■ ■ • x X n must be based, to start with, on a formulation on each factor X y 
The Hardy space on Xj is then defined by developing the Littlewood-Paley theory on X y 
Our approach includes the following five steps. 

1. Introduce the spaces of test functions and distributions. I 11 the classical case, 
the relevant spaces of test functions and distributions are just Schwartz test functions and 
the class of tempered distributions. I 11 order to study the Calderon reproducing formula 
associated with the T(b) theorem, the new test function and distribution spaces were first 
introduced on Euclidean spaces in (HU, and on spaces of homogeneous type, where the 
quasi-metric d satisfies the Holder regularity condition (jl.4[) and the measure /i satisfies 
condition (11.51) . in [ HS] . See [HMY J and I I 11.1.2 . respectively, for spaces of test functions and 
distributions on spaces of homogeneous type with additional assumptions. In this paper, we 
will introduce test functions and distributions on spaces of homogeneous type in the sense of 
Coifman and Weiss. These spaces include all those considered previously. 

2. Establish the wavelet reproducing formula on test functions and on dis¬ 
tributions. The classical Calderon reproducing formula was first used by Calderon in 0- 
Such a reproducing formula is a powerful tool, particularly in the theory of wavelet analy¬ 
sis. See (MU. Using Coifman’s decomposition of the identity operator, as mentioned above, 
David, Journe and Semmes [ DJS] gave a Calderon-type reproducing formula which was a key 
tool in proving the T(b) theorem on M n and the T(l) theorem on spaces of homogeneous type 
with the conditions (jl .4[) and (11.51) . See [ HMY J and [HLL2] for the continuous and discrete 
Calderon reproducing formulas on spaces of homogeneous type with additional assumptions. 
As mentioned above, Auscher and Hytonen established a wavelet expansion on L 2 ( X) (and on 
L P (X), 1 < p < 00 ). For our purposes, we will show that the wavelet expansion constructed 
in EU also converges in both the test function and distribution spaces. 

As Meyer pointed out in |Mlj . “The wavelet bases are universally applicable: ‘everything 
that comes to hand’, whether function or distribution, is the sum of a wavelet series and, 
contrary to what happens with Fourier series, the coefficients of the wavelet series translate 
the properties of the function or distribution simply, precisely and faithfully. ” I 11 particular, 
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our results provide such wavelet expansions for test functions and for distributions, and are 
used below to introduce square functions and develop the Littlewood-Paley theory. 

3. Develop the Littlewood—Paley theory. Based on the wavelet expansion provided 
in era. one can formally introduce two kinds of square functions, namely, the discrete 
version defined in terms of wavelet coefficients and the continuous version defined via wavelet 
operators D & (different from the operators Dj. = Sk+i — Sk mentioned above). To show that 
the L p norms of these square functions are equivalent, for a suitable range of p, we need a 
Plancherel-Polya inequality. The classical Plancherel-Polya inequality says that the L p norm 
of a function / whose Fourier transform has compact support is equivalent to the £ p norm 
of the restriction of / to an appropriate lattice. This kind of inequality was first proved in 
jH2j on spaces of homogeneous type with the conditions (11.41) and (11.51) . and in jHMYj and 
[HLL2], respectively, for the one-parameter and multiparameter cases with some additional 
assumptions. As a consequence of our Plancherel-Polya type inequalities, the Hardy space on 
spaces of homogeneous type in the sense of Coifman and Weiss is well defined. In particular, 
as in the classical case, H p = L p for 1 < p < oo. 

4. Introduce the generalized Carleson measure space. It is well known that in 
the classical one-parameter case, the space BMO, as the dual of II 1 , can be characterized by 
Carleson measures. Moreover, in [CFl] Chang and Fefferman proved that the dual of product 
H 1 is characterized by product Carleson measures. The generalized Carleson measure space 
CMO p , as the dual of the product H p , was introduced in | iHLL1i | and [ HLL2 ] on spaces of 
homogeneous type with some additional assumptions. In the current paper, working in the 
setting of spaces of homogeneous type in the sense of Coifman and Weiss with no additional 
assumptions, we introduce CMO p in terms of wavelet coefficients, and prove that CMO p is 
the dual of H p . In particular, CMO 1 = BMO is the dual of H l . Moreover, we also introduce 
the space VMO and show that VMO is the predual of H l . 

5. Establish the Calderon—Zygmund decomposition. The classical Calderon- 
Zygmund decomposition played a crucial role in developing Calderon-Zygmund operator 
theory. This decomposition has many applications in harmonic analysis and partial differ¬ 
ential equations. Such a decomposition for product Euclidean spaces was first provided by 
Chang and Fefferman in (CF2j . The main tool used in [CF2] is the atomic decomposition. 
In the current paper, applying the wavelet expansion constructed in [ AH] , we establish the 
Calderon-Zygmund decomposition on product H p on spaces of homogeneous type with no 
additional assumptions. As a consequence, we obtain the interpolation of operators that are 
bounded from Hardy spaces to Lebesgue spaces, and of operators that are bounded on Hardy 
spaces. 

We note that in the original work on extending the Calderon-Zygmund theory to spaces of 
homogeneous type (A", d', p), the philosophy was as follows: Coifman constructed the approx¬ 
imations to the identity Sk, which were used in [D.TSj to define the continuous square function 
and to establish the Littlewood-Paley theory. Later the discrete Calderon reproducing for¬ 
mula was introduced and the Littlewood-Paley theory for the classical function spaces were 
established in ED and IHSI, respectively. By contrast, in our setting of (X,d,p) with the 
original quasi-metric d, we begin with the discrete wavelet reproducing formula (Theorem 13.4]) 
and define the discrete square function S(f) in terms of wavelet coefficients (Definition 14. ip . 
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However, there does not seem to be a direct proof of the Littlewood-Paley theory for S(f). 
The question then is: how to find a continuous version of the square function? We introduce 
a new continuous square function S c (f) (Definition 14.21) . via certain wavelet operators Z\. 
that are expressed in terms of the Auscher-Hytonen wavelets (Lemma I3.6[) . We prove that 
||S'c(/)|| P ~ ll/llp for 1 < p < oo (Theorem [OJ, and that ||S(/)|| P ~ ||S' C (/)|| P both for 
1 < p < oo and moreover for an additional range of p < 1 depending on the upper dimen¬ 
sions lO i of the factor spaces X ^ and on the Holder regularity exponents ry of the wavelets 
(Theorem 14.3j) . 

We remark that in this paper we concentrate on the product case. As Nagel and Stein 
observed in (NSlj . “Any product theory tends to be burdened with notational complexities.” 
For notational simplicity, we have written our results and proofs for the case of two param¬ 
eters. However, our methods also establish the corresponding results for the product case 
with k factors, for fceE We also point out that these results extend related previous results 
from |D.TS1 IH11 IH2[ IHLL1[ IHLL2[ IHMY1 IHS] and the references therein. In those papers 
either extra assumptions are made on the quasi-metric and the measure, or the product case 
is not covered, or both. 

The paper is organized as follows. In Section [2] we briefly recall the systems of dyadic cubes 
from [ HK ] and the orthonormal bases from [ AH ] on spaces of homogeneous type in the sense of 
Coifman and Weiss. In Section [3] we introduce the one-parameter and product test functions 
in Definitions 13.11 and 13.91 respectively, together with the induced classes of distributions. 
The main result in this section is Theorem 13.41 which gives the wavelet reproducing formula 
for test functions. In Section UJ the Littlewood-Paley square functions in terms of the wavelet 
coefficients and of the wavelet operators are given in Definitions 14.11 and 14.21 respectively. 
The two main results here are Theorems 14.31 and 14.41 Theorem 14.31 gives the Littlewood- 
Paley theory, including the norm equivalence of the discrete and continuous Littlewood-Paley 
square functions. Theorem 14.41 gives the Plancherel-Polya inequalities, which are the main 
tool in proving Theorem 14.31 The product H p , CMO p , BMO and VMO spaces are defined 
in Section [5] via the orthonormal wavelet basis. We use the Plancherel-Polya inequalities 
again to show that these function spaces are well defined. The duality results are given in 
Theorem 15.31 for ( H p )' = CMO p and in Theorem 15. 101 for (VMO)' = H 1 . Finally, in Section[6] 
we prove the Calderon-Zygmund decomposition and the interpolation theorem for Hardy 
spaces in Theorems 16.11 and 16.21 respectively. 

2. Preliminaries 

We are interested in establishing the Hardy space theory on spaces X — X x x X 2 x • • • x X n . 
Each factor is a space of homogeneous type in the sense of Coifman and Weiss. We will first 
need to develop a Littlewood-Paley theory for each factor X i} 1 < i < n, and then pass to 
the corresponding product theory. In this paper, we always assume that p(Xi) = oo and 
that fj,(B(x,r)) > 0 for all r > 0 and x 6 X*, for 1 < i < n. As usual, C denotes a constant 
that is independent of the essential variables, and that may differ from line to line. 

In this section we recall the systems of dyadic cubes, in a geometrically doubling metric 
space, constructed by Hytonen and Kairema [HK]; and the orthonormal wavelet basis, on 
spaces of homogeneous tvpe. constructed bv Auscher and Hvtonen (AH l SAH2]. See also EE). 
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[ AH'!] and the references therein for the history and applications of various versions of dyadic 
cubes. 

2.1. Systems of dyadic cubes in a geometrically doubling metric space. Let d be 

a quasi-metric (defined in the Introduction) on a set A". The quasi-metric space (A, d) 
is assumed to have the following geometric doubling property: there exists a positive integer 
Ai E N such that for each x E X and each r > 0, the ball B(x, r) := {y E X : d(y, x) < r} can 
be covered by at most Ai balls B(xi,r/ 2). It is shown in [CW1] that spaces of homogeneous 
type have the geometric doubling property. 

As usual, a set h C A" is open if for every i6fl there exists e > 0 such that B(x, e) C fi, 
and a set is closed if its complement is open. 

Theorem 2.1 ([HK] Theorem 2.2). Suppose that constants 0 < Co < Co < oo and 6 E (0,1) 
satisfy 

( 2 . 1 ) l2AlC 0 5 < co. 

Given a set of points {z k } a , a E sAk, for every k E Z, with the properties that 

(2.2) d(z k z k g ) > c 0 S k ( a ^ (3), min d(x,z k ) < C 0 5 k , for all x E X, 

^ a 

— , . _/j. 

we can construct families of sets Q k a C Q k a C Q a (called open, half-open and closed dyadic 
cubes,), such that: 

(2.3) Q k a and Qjf are the interior and closure of Q k a , respectively, 

(2.4) if £ > k, then either Qp C Q k a or Q k a n Qp = 0; 

(2.5) A = ( disjoint union ) for all k E Z; 

a 

(2.6) B(z , f,c 1 5 k ) CQ‘C B(z k , Ci5 k ), where c 1 := (3Ao) _1 c 0 and C x := 2A 0 C 0 ; 

(2.7) if l > k and Q J C Q k a , then B^Crf*) C B^C^). 

The open and closed cubes Q’f and Q k a depend only on the points Zp for t > k. The half-open 
cubes Q k a depend on Zp for £ > min(A:, k 0 ), where k 0 E Z is a preassigned number entering 
the construction. 

2.2. Orthonormal wavelet basis and wavelet expansion. In this subsection, we recall 

the orthonormal basis and wavelet expansion in L 2 (A) which were recently constructed by 
Auscher and Hytonen [AH]. To state their result, we must first recall the set { x of reference 
dyadic points as follows. Let 5 be a fixed small positive parameter (for example, as noted 
in Section 2.2 of (AH], it suffices to take 6 < 10 -3 Ao 10 ). For k = 0, let := be a 

maximal collection of 1-separated points in A". Inductively, for k E Z + , let & k := { xD 

k ~ l and %~~ k := C p e § k _ anc [ 5 _fc -separated collections in and 

respectively. 

Lemma 2.1 in [ AH] shows that, for all k E Z and x E X, the reference dyadic points satisfy 

(2.8) d(x k a ,x k p) > 5 k (a ^ (3 ), d(x, 3T k ) = min d{x,x k a ) < 2A 0 5 k . 

^ a 
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Also, taking Co := 1, Co := 2A 0 and 5 < 10 ” 3 Aq 10 , we see that Co, Co and 5 satisfy (12. ip in 
Theorem 12.11 Therefore we may apply Hytonen and Kairema’s construction (Theorem 12.ip . 
with the reference dyadic points {Xa}k&z,a£^ k playing the role of the points {z k }kez,ae£ 4 , to 
conclude that there exists a set of half-open dyadic cubes 

{Qa}kez,ae$r k 

associated with the reference dyadic points {Xa}k£Z,ae%' k - We call the reference dyadic 
point x k the center of the dyadic cube Q k . We also identify with the set of indices a 
corresponding to x k a G 

Note that C £C k+l for fc G Z, so that every x ^ is also a point of the form x k C~ l . We 
denote ( 3L k := £C k+1 \^C k , and relabel the points {x^} a that belong to & k as {Va}a- 

We now recall the orthonormal wavelet basis of L 2 (X) constructed by Auscher and Hytonen. 


Theorem 2.2 (|AH] Theorem 7.1). Let (X,d,p) be a space of homogeneous type with quasi¬ 
triangle constant A 0 , and let 

(2.9) a:=(l + 21og 2 A 0 )- 1 . 


There exists an orthonormal wavelet basis k G Z, y k G ( 3L k , of L 2 (X), having exponen¬ 

tial decay 


( 2 . 10 ) 


W(*)l < 


C 


Vh(B(y k ,8 k )) 


exp 


d(y k ,x)y 


5 k 


Holder regularity 

(2.11) |^(*)-^(y)|< 


C 

Vh(B(y k ,8 k )) 


d(x,y)y / ft 

eXP (-K- 


(d(y k ,x)y 


8 k 


for d(x, y) < and the cancellation property 


( 2 . 12 ) 



v£(z) d h(x) = 0, 


for ke Z, y k a G W k . 


Moreover, the wavelet expansion is given by 

(2.i3) = 

k£Z 

in the sense of L 2 (A"). 

Here 5 is a fixed small parameter, say S < 10~ 3 Aq 10 , and C < 00 , v > 0 and rj G (0,1] are 
constants independent of k, a, x and y k . 

In what follows, we refer to the functions f) k x as wavelets. Throughout this paper, a denotes 
the exponent from (12.91) and r) denotes the Holder-regularity exponent from (12.111) . 


Remark 2.3. The wavelets {f> k }k,a form an unconditional basis of L P (X) for 1 < p < 00 , 
as shown in Corollary 10.4 in [AH]. Therefore, the reproducing formula (12.131) also holds 
for / G L P (X). Moreover, for us the most crucial feature of the orthonormal wavelets 
construction is the following estimate, which is a special case of |AH1 Lemma 8.3]: 


(2.14) 


£ 

j€Z:Si>r 


1 

li(B(x,6i)) 


exp 


/ d(x,& j ) y\ C 

V 8? ) ) ~ /a(B(x,r))' 
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for all x G X and r > 0, and for the constants v > 0 and a := (1 + 21og 2 Ao) 1 from 


Theorem E2J Series of this type naturally arise in the context of proving that the reproducing 
formula holds for test functions and distributions, as well as in relation to function spaces. 
This estimate allows us to drop the extra assumption, used in previous work, of a reverse- 
doubling property on the measure. 

Furthermore, this estimate is crucial for estimating the quantity ^a( x )^a(y)- We 

will use this estimate in the proofs of Theorems 13.41 and 14.31 below. 

3. Test functions, distributions, and wavelet reproducing formula 

We now introduce test functions and distributions on spaces of homogeneous type (. X, d, //) 
in the sense of Coifman and Weiss, and on product spaces (Ad, d\, pi) x (X 2 , d 2 , y 2 ). We show 
that the (scaled) Auscher-Hytonen wavelets are test functions (Theorem 13.3p . and establish 
the wavelet reproducing formula for test functions and for distributions (Theorem 13.41 Corol¬ 
lary 13.51 Theorem 13.1 111 . Along the way we establish some properties of the wavelet opera¬ 
tors Dk (Lemma 13.61) . and construct smooth cut-off functions using the splines of Auscher 
and Hytonen (Lemma 13.81) . 

We begin with the one-parameter case. For x, y G X and r > 0, let 

V r (x) n(B(x, r)) and V(x,y) := y(B(x,d(x,y))). 

3.1. One-parameter test functions, distributions, and wavelet reproducing for¬ 
mula. 

Definition 3.1. (Test functions) Fix xq G A", r > 0, f3 G ( 0 , 77 ] where 77 is the regularity 
exponent from Theorem 12.21 and 7 > 0. A function / defined on A" is said to be a test 
function of type (x 0 ,r,/3, 7 ) centered at x 0 G A if / satisfies the following three conditions, 

(i) (Size condition) For all x G X, 



(ii) (Holder regularity condition) For all x, y G X with d(x,y) < (2A 0 ) 1 (r + d(x,x 0 )), 



(iii) (Cancellation condition) 

f(x) dn(x) = 0 . 



A priori , this definition makes sense for arbitrary /3 > 0. Here we have used the condition 
/3 G (0, 7 ] both for consistency with the earlier literature and since our focus is on the 
wavelets if*, which (when scaled) are test functions with /3 = 77, as we will see. 

We denote by G(x 0 , r, /3, 7 ) the set of all test functions of type (x 0 , r, ft, 7 ). The norm of / 
in G(xo,r,j3, 7 ) is defined by 


II/IIg^oxAt) := inf {C > 0 : (i) and (ii) hold}. 


For each fixed xq, let G(/3, 7 ) := G(xq, 1, ft, 7 ). It is easy to check that for each fixed 
x'q G X and r > 0, we have G{x' Q , r, ft, 7 ) = G(J3, 7 ) with equivalent norms. Furthermore, it 
is also easy to see that G(j3, 7 ) is a Banach space with respect to the norm on G(ft, 7 ). 
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For (3 £ ( 0 , ■//] and 7 > 0 , let G(P, 7 ) be the completion of the space G(rj, 7 ) in the norm of 

o o 

G(f3, 7 ); of course when (3 = r/ we simply have G(P, 7 ) = G{j],t) = G(r), 7 ). We define the 
norm on G(P, 7 ) by 11/11°^ : = II/IIg(/3, 7 )- 

O 

It is immediate from the definition that the sets G(P, 7 ) are nested; for example if 0 < 
P' < (3 and 0 < 7 ' < 7 , then G{P,l) C G{P', 70- 


Definition 3.2. (Distributions) Fix xq £ X, r > 0, P £ ( 0 ,77] where rj is the regularity 

O 

exponent from Theorem 12.21 and 7 > 0. The distribution space (G{P, 7))' is defined to be 

O 

the set of all linear functionals C from G{P, 7) to C with the property that there exists C > 0 

O 

such that for all f £ G{P, 7), 

IU/)I< 01/11. (ft7) . 

We note that for each p £ ( 0 , rj] and 7 > 0 , the set G(P, 7) C L 2 (X ), while each / £ L 2 (X ) 

o 

induces a distribution in (G(P, 7))'. 

We now prove that the wavelets constructed in [XH], suitably scaled, are test functions. 

Theorem 3.3. Suppose that {V’a}*:ez,ae® rfc an orthonormal wavelet basis as in Theo- 
rem 12.21 with Holder regularity of order rj. Then for each k £ 7L, y k £ ( 3V k , and 7 > 0, 
the scaled wavelet ip k (x) / \Jfi(B(y k , S k )) belongs to the set G(y k ,5 k ,rj,'f) of test functions of 
type (y k , 5 k , rj, 7 ) centered at y k . 


Before proving Theorem 13.31 we make the following useful observation: by the doubling 
property (jL3j) on the measure p, for each x 0 , x £ X and r > 0 with r < d(x 0 ,x), we have 
V(xq,x) < C(d(x 0 ,x)/r) ul V r (x 0 ) and hence 


(3.1) 


1 < c 1 _ r dfx {h x) y 

V r (x 0 ) ~ V r (xo) + V{xq,x) \ r ) 


Proof of Theorem Iff. <51 By property (12.10p of r f k from Theorem 12.21 we obtain that 

4>a( x ) - C 


Vp( B (Va’ Sk )) 


< 


< 


exp(-z ;(5- k d(y k ,x)) 


ji(B(y k ,5 k )) 

C / S k 


V S k(y k ) + d(y k ,x), 

for all T > 0, with a constant C depending only on u, a = (1 + 21og 2 ffl 0 )~ 1 , and T. To see 
that ifa( x )/\/h(B(y k , 5 k )) satisfies the size condition Definition Mi) for 7 > 0 , we consider 
two cases. First, if 5 k > d(y k ,x), then V(y k ,x ) < V$k(y k ) and hence 

i’afa) < _ G _/ 

^/ji(B(y k ,5 k )) ~ V 5 k(y k ) + V(y k ,x)\S k + d(y k ,x) 

For the second case, if 5 k < d(y k ,x), an application of (13. ip shows that 




< C- 


r-w 


yjji(B(y k ,5 k )) V S k(y k ) + V(y k ,x)\6 k + d(y k ,x 

Taking T > ui and setting 7 := T — 00 , we see that f} k a {x)/yjp(B(y k , 8 k )) satisfies the si 
condition Definition 13.lT i) with xq = y k and r = 5 k , and for arbitrary 7 > 0. 


size 
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We now show that A( x ) / \Jy(B(y k , S k )) satishes the smoothness condition Defiriition l3.in i) 
for 7 > 0 and j3 = rj, if d(x,y ) < (2A 0 )^ 1 (5 k + d(y k ,x)). We consider three cases. First, 
suppose d(x,y) < 5 k < (2 A 0 )~ 1 (S k + d(y k ,x)). Then property (12.111) yields 


A a ( x ) ~ i>a{y) 

< C 

( d{x,y)\ 

v 1 

( ^ ) 

VA B (y k ,5 k )) 

[ 5 k > 

»(B(y k ,6 k )) 

\8 k + d(y k , x) / 


Note that in this case, 5 k < (2A 0 — 1) 1 d(y k , x), and so we may apply (13.11) with r = 5 k and 
xq = y k to conclude that 


a ( X ) - ^a(y) 

^ c ( d ( x iV) \ 

v 1 / 5 k \ 

VA B (y k ,S k )) 

\S k + d(y k ,x)J 

v s^(ya) + y(yt x ) \8 k + d(y k ,x)j 


Second, consider the case where S k < d(x,y ) < ( 2A 0 ) 1 (8 k + d(y k , x)). It is straightforward 
to verify from the quasi-triangle inequality that in this case, 

(3.2) 5 k + d(y k ,y)~5 k + d(y k ,x). 


Applying property (12.101) for Aix ) and ^(y), we hnd 


A a ( x ) ~^a(y) 

<r C 

/ 8 k 

Y+ c 

( 6k Y 

Vv(B(y k ,S k )) 

1 / 

e s- 

+ d(y k , x) 

) vAy k ) 

\ 5 k + d(y k ,y)J 


< c ( d{x,y) y 1 / 5 k \ r 

\ 6 k J V S k(y k )\6 k + d(y k ,x)J 

< / d{x,y) y _ 1 _/ $ k y-u-v 

\S k + d(y k ,x)J V 5 k(y k )+ V(y k ,x)\5 k + d(y k ,x)J 

Here the second inequality follows from (13.21) and the fact that d(x,y)/S k > 1, and the third 
inequality follows from (13. ip . 

For the third and last case, if S k > (2A 0 ) _1 (h fc + d(y k , x)) > d(x, y), then we have d(x, y) < 
5 k and d(y k ,x) < (2 A 0 — 1)<W Therefore, applying property (12.lip together with the fact 
that V(y k ,x) < Cy(B(y k ,8 k )) yields 


Tpa( X ) ~ $a(y) 

< c 

(d{x,y)\ 

” 1 

f s " \ 

VA B (y k ,8 k )) 

{ 5 k J 

vAvi) 

Kd k + d(y k ,x)J 


< c ( d ( x iv) y_i_/ $ k 

\ 5 k + d(y k ,x)J V S k(y k )+ V(y k ,x)\ 6 k + d(y k ,x) 


Combining all the cases above, in the hrst and second cases take T > uj + 77 and set 
7 := T — u> — y, and in the third case take V > y and set 7 := T — rj. We see that the function 
k (x )/ \Jy(B(y k , S k )) satishes the smoothness condition (ii) in Definition 13.11 with x 0 = y k , 
r = 6 k , /3 = r) and for arbitrary 7 > 0. 

Moreover, the cancellation property of the function ip k (x)/y/y(B(y k , S k )) is immediate 
from property (12.121) of ij} k , by Theorem 12.21 Thus, ip k {x)/\Jy(B(y k , 5 k )) belongs to the test 
function space G(y k , 5 k , 7 , 7 ). This completes the proof of Theorem 13.31 □ 


Now we state and prove the main result of this subsection, which will be the crucial tool 
for establishing the Littlewood-Paley theory and developing the Hardy spaces. 
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Theorem 3.4. (Wavelet reproducing formula for test functions) Suppose that f G Giftft 
with ft 7 G (0,?/). Then the wavelet reproducing formula 

(3.3) /w = EE</-*mw 

fcez ae^ k 

holds in G(ft, ft) for all ft G (0,/3) and ft G ( 0 , 7 ). 

As an immediate consequence of Theorem 13.41 we obtain the following result. 


Corollary 3.5. (Wavelet reproducing formula for distributions) Take ^76 (0 ,rj). Then 

O 

the wavelet reproducing formula (j3.3[i also holds in the space (Gift ft)' of distributions. 


Proof of Theorem\3.f\ Take / G Gift ft with ft 7 G ( 0 , 7 /). It suffices to show that 


(3.4) 


C(D',V) 


0 


E E^‘-/)v>7) 

|fc|>L 

as L tends to infinity, for each ft G (0, ft and ft G ( 0 , 7 ). 

The proof of (13.41) is based on the following estimate: for each ft G (0, /3) and ft G ( 0 , 7 ), 
there is a constant cr > 0 such that for each L G N, 


(3.5) 


E E «./«<•) < 


where C is a positive constant independent of / G Gift ft- 

To verify (13.5p . it suffices to show that the following decay and smoothness estimates hold: 
for each ft G (0, ft, there exist a positive constant C independent of /, and a positive number 
cr', such that 

1 / 1 \i 

|g(/3 i7 ), and 


( 3 . 6 ) | e E M 

\k\>LaeW k 


[X 


< Cd a ' L 


Vftxo) + V(x,x 0 ) VI + d(x,x 0 ) 


( 3 . 7 ) | e E -EE 

\k\>La£ty k \k\>La&W k 

/ rffas') W 1 / l y ||n| 

— Vl + d(x, x 0 ) / hi(^o) + V(x,Xq) Vi + d(x, x 0 )J G ^,i) 
for all x and x' such that d(x,x') < ( 2 A 0 ) _1 (1 + d(x,x 0 )). 

Indeed, to see that (13.6[) and (j3.7[) imply (13.5p . we take the geometric mean between (13. 7[) 
and the following estimate (obtained directly from 


E E»t/)*)-E E«£v)^m 


|fc|>L a£$f k 


< 


\k\>La& 3 f k 


< cs 


E E «£• m o 

\k\>L 

1 


[X 


cr' L 


Vftxftj + V(x,xft Vl + d(x,x 0 ) 


E E«’4/)^m 

\k\>L a^ 9 k 

1 ^ 7 ' 

" "‘1008,7) 


This gives 
(3.8) 


E E w tMio) -EE 
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< CS 


<jL 


d(x, x') 


<r + d(x,x o) 

for some a < a'. Now (13.61) and 


Vi (xo) + V(x, x Q ) V 1 + d(x, x Q ) 
, together with the fact that 


II/IIg , (/3,7) 



E (V£,/}V£(®) d l*(x) = °> 


imply that '52 a e& k (' t l } an f)^a is a test function in G(/3', 7 '). Moreover, we see from the 

upper bounds in (13.61) and (13.81) that (13.51) holds, as required. 

To prove the decay and smoothness estimates (13.61) and (13.7f) . we need the following lemma 
which gives estimates for the kernels D k (x,y ) = ipa( x )ipa{y) of the wavelet opera¬ 

tors D k . These wavelet operators are defined by 

D k f(x) := f)^a( x ) = f D k (x, y)f(y) dy, 

Q.ggl'fc ^ 

for k G Z. We note that the first two estimates in Lemma 13.61 are similar to estimates given 
in Lemma 9.1 in [AH ]. 


Lemma 3.6. (Properties of wavelet operators D k ) Let 

D k (x,y) := E 

a£& k 

for x, y e X. Fix 7 > 0. Then the following estimates hold. 
(i) (Decay condition) For all x, y G X, we have 

1 / 5 k 


(3.9) 


\D k (x,y)\ < C- 


V S k (x) + V(x, y) V S k + d(x, y) < 
ii) (Smoothness condition) Ifd(y,y') < ( 2 A 0 ) _1 max{<5 fc + d(x, y), S k + d(x, y')}, then 


(3.10) \D k (x,y) - D k (x,y')\ < c(- 


d{y,y') 


(3.11) 


S k + d(x,y)J V 5 k(x) + V(x,y) \5 k + d(x,y ) 

The same estimate holds with x and y interchanged. 
iii) (Double smoothness condition) If 

d(x, x') < (2A 0 ) _1 max{h fc + d(x, y),5 k + d(x ', y)} and 
d(y , y') < ( 2 A 0 ) _1 max{h fc + d(x, y), 5 k + d(x, y')}, 

then 

| D k (x,y) - D k (x',y ) - D k (x,y') + D k (x',y')\ 

< C 


\V d (y,y') ) 

i" 1 1 

f i 

' \5 k + d(x,y )/ 

V S k(x ) + V(x,y) 

V 5 k + d(x, y)J 


We defer the proof of Lemma 13.61 until after the end of the proof of Theorem 13.41 
Returning to the proof of Theorem 13.41 we first show (13.61) . Write 


E E<^-/></>; 

|fc|>L 


[X 


< 


E E M’M 

k>L oiG^k 

= : (A) + (B). 


lx 


E E^>^) 

k<—L a£ <3/ k 
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For (A), using the cancellation property (12.12j) of the wavelet yj k and integrating over the 
sets W\ := {y E X : d(x,y) < (2A 0 ) _1 (1 + d(x,x o))} and W 2 : = X \ W\ , we obtain 

( A ) < ^2 f \Dk(x,y)\ \f(y) - f(x)\ dfi(y) 

k>L JW 1 

+ 5^ [ \D k {x,y)\{\f(y)\ + \f(x)\)dy(y) 

k>L JW * 

= '■ ( A )l + ( A )2- 


To deal with (A) 1 , applying the decay condition (13. 9 p from Lemma [3761 on D k (x,y ) and the 
Holder regularity property (Definition I3.1f ii)) of the test function / G G(/3, 7 ) gives 

(A), < Vst ( x) + V(i,!/) G* + 

/ d(x,t/) x/ 3 1 / 1 XT' 

VI + d(x, x 0 )J Vi(x 0 ) J rV(x,Xo)\l +d(x,x 0 )J 

~ ^ Vi(xo) + V(x,x 0 ) (l + d(x, x 0 )) 

To estimate (A) 9 , applying the size conditions on both D k (x, y) (Lemma l3.6p and / (Dehnition 
Mi)) gives 


( a ) 2 < c 'll/llG(/3, 7 )^y^ v sk ( x ) + V(x,y) \5 k + d(x,y) 


k>L 

.V'i(xo) + V(y,x 0 ) (l + d(y,x 0 )) + V^Xq) + V(x, ar 0 ) (l + d{x, x 0 )) 

For the hrst sum, use the fact that if d(x,y) > (2A 0 ) _1 (1 + d(x,x 0 )) then V(x,y) > 
CV(x, (2A 0 ) _1 (1 + d(x, x 0 )) > C\Vi(xq) + V(x, x 0 )]. For the second sum, apply the following 
estimate: 


We obtain 


'w 2 v 6*(x) + V(x,y)\5 k + d(x,y) 

1 


M < cs^ k . 


(A) 2 <C\\f\\ G(0 ^ L 


Vx{x Q ) + V(x,x 0 ) VI + d(x,x 0 ) 


Hence for some a > 0, 

(A) < C6 aL 


Vi(x 0 ) + V(x,x Q ) VI + d(x,x 0 ) 


II/IIg(0,7)- 


We now turn to (B). Using the fact that f x f(x)dy(x) = 0 and considering the sets 
W 3 := {y E X : d{y,x 0 ) < (2A 0 )~ 1 (5 k + d{x,x 0 ))} and W 4 := X \ W 3 , we have 

(B) < f \ D k( x >y) -Dk(x,x 0 )\\f{y)\dfi{y) 

k<-L Jw 3 

+ yz l (\ D k(x,y)\ + \D k {x,x 0 )\)\f(y)\dy(y) 

k<-L Jw ± V ' 

= : (B) 1 + (B) 2 . 
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For (B) 1 , applying the smoothness estimate from Lemma [3.6I ii) and the size estimate of 
the test function / (Definition 13. If i)) yields 


(B)l £ Gl/lleUW £ Jw 3 (s k + d(x!x 0 )) V s ,(x) + V(x,x „) 


X 


k<-L 
6 k 


5 k + d(x, x 0 )) V\ (x 0 ) + V(y, x 0 ) V1 + d(y, x 0 ) 

1 /I 


dy(y) 


< CWlWa^^, 


Vl (x 0 ) +V(x,x 0 )\1 + d(x, x 0 ) - 
where we choose 7 1 < rj' < 7 . To estimate (B) 2 , we first write (B ) 2 = (B ) 21 + (B ) 22 where 


(B) 21 : = E / \ D k(x,y)\\f(y)\dfi(y) 

k<-L Jw ± 


and 


( B ) 22 := X] / \ D k{x,x 0 )\\f(y)\dn(y). 

k<-L Jw ± 


Since here d(y,x 0 ) > (2A 0 ) 1 (<5 fc + d(x,x 0 )), the size estimates for the test function / (Defi¬ 
nition I3.1f i)) imply that for 0 < 7 ' < 7 , 

'/<*')' £ C ' II/IIg1 ' , ' 7) n(zo) + v(y, x 0 ) (1 + d(y, x a ) 

< CII/llGowi*^) 1 


Vl (x 0 ) + V(x, x 0 ) V1 + d(x, x 0 ) ■ 

The above estimate, together with the fact that f Y \D k (x,y) \ dy(y) < C, yields 

1 


(B) 21 < C8^' )l , 


II/IIg(/ 3 , 7 )- 


Vi(x 0 ) + V{x,x Q ) Vl + d(x,x 0 ). 

The estimate for (B ) 22 is similar, but easier. Indeed, by Lemma [3.61 we have 


Ini \ \ ^ r* 1 ( 8 k ^ 7 ' 

\D k x,x 0 | _ ^yk^ XQ ^ _)_ v(x,x 0 ) V(5 fc + d(x,x 0 )/ 

and 

Thus, we obtain the same estimate for (B ) 22 as for (B) 21 , but with 7 — 7 ' replaced by rf — 7 '. 
This completes the proof of (13.61) . 

Finally, we show (13.71) . To do this, we first need to construct a smooth cut-off function. 
For this purpose, we recall the following result on the properties of the spline functions s 
on (X,d,fi) that were constructed by Auscher and Hytonen [ AH] . 


Theorem 3.7 ([AH], Theorem 3.1). The spline functions s * satisfy the following properties: 
bounded support 

XB(x*,l/8A~ 3 S k )( X ) — S a( X ) — XB(a;fc,8A^(5 fc ) (x), 
the interpolation and reproducing properties 

4(xj) = <w, = 1 > s « = 

a /} 


P 
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where {p^g}p a finite nonzero set of nonnegative coefficients with p k a g < 1; and Holder 
continuity 


s k a {x) - s k a (y)\ < C 


( d(x,y) y 


We point out that in the above theorem, a runs over 3C k . Using these splines we can 
construct a smooth cut-off function, as follows. 


Lemma 3.8. (Smooth cut-off function) For each fixed x 0 £ A" and R 0 £ (0, oo), there exists 
a smooth cut-off function h(x) such that 0 < h(x) < 1 , 

h(x) = 0 when x £ B(x 0 , Agi? 0 ) c , 


d(x,y)y 
5 k ) ' 


h(x) = 1 when x £ B(x 0 , R 0 /A), 

and there exists a positive constant C independent of x o , Rq, x, y such that 

\h(x)-h(y)\<C 

Proof. For a hxed Ro £ (0, oo), we choose ko £ Z such that 

8Ag(5 fc ° < Rq/A and > R 0 / 4. 

Next, we dehne the index set X k(i as follows: 

l ko := {a £ : B(x%, 8A 5 0 S k °) n B(x 0 , Ro/4) ^ 0}. 

Then the number of indices contained in X ko is bounded by a constant independent of R 0l 
k 0 , and x 0 , since 8Ag<5 fc ° is comparable to R 0 and the reference dyadic points {ic^, 0 } are 
<5 fc °-separated. 

Now dehne 

h(x) := s *( x )- 

OL^fil-kQ 

From the properties of the spline functions s k (x) (Theorem 13.7jh it is easy to verify that h(x) 
satishes all the properties listed in Lemma [3.81 □ 


In what follows, the Holder-regularity index of the cut-off function h(x) is the 77 given in 
Theorem 13.71 ([AH]. Theorem 3.1). 

We now return to the proof of (13.71) . It suffices to show that there exists a constant C such 
that for each M > 0 and for d(x,x') < (2H 0 ) _1 (l + d(x,x 0 )), 


(3.12) 



[D k (x,y) - D k (x',y)\f(y) dp(y) 


/ d(x,x') \ 

f 1 ( 1 ) 

Vl + d(x,x q)J 

Vfixo) + V(x,xf) Vl + d(x,xo)J 


To do this, fixM >0 and let T denote the wavelet operator given by 


T (/) 0 ) : = [ K{x,y)f{y)dn{y), 

Jx 

with kernel K(x, y) := D k (x, y), where the kernel D k (x, y) of the wavelet operator D k 

is given by ■ 
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To show (I3.12j) . it suffices to prove that if d(x,x') < (2A 0 ) _1 (l + d(x,x o)), then 

|r ( / )W - T(/ )M | < c( T ^^) V lW+ V (j , Ja) ( TT ^y) 7 - 

Let R = d(x, xq) and r = d(x, x'), and consider the case where R > 10 and r < (20Ao) _1 (l + 
d(x,x o)). Following |M2j . set 1 — I(y) + J(y ) +L(y), where I(y) is a smooth cut-off function 
as in Lemma [3.81 satisfying 

I(y) = 1 when y G B(x, R/32AI) and /(?/)= 0 when y 6 B(x, R/8) c , 

and 

J(y) = 1 when y G B(x 0 , R/32A 2 0 ) and I(y) = 0 when y e B(x 0 , R/8) c . 

Also set 

/l(y) := f(y)I{y), f*{y) ■= f(y)J(y ) and f 3 (y) := f(y)L(y). 

It is easy to see that /i, / 2 and f 3 satisfy the following estimates fid. 1.31) (13.171) : 

1 / 1 


(3.13) 


\fi(y)\ < C-||/|| goj>7) : 


Vi (x 0 ) + V(x, Xq) V 1 + d(x, X 0 ) 
since \fi{y)\ < \f(y)\ and 1 + d(y,x 0 ) > C(1 + d(x,x 0 )) by the form of /i; 

\h(y) - h(y')\ < \imf{y)-fW)\ + \fW)my)-iW)\ 

(3-14) < C\\f\\ Gi M 

for all y and ?/; 

(3.15) |/ 3 (y)| <011/110^) 


f d(y, y') \ 

\ 1 ( 1 ) 

Vl + d(x, x 0 )) 

Vi(x 0 ) + V(x,xq) Vl + d(x,x 0 )J 


(3.16) 
and 

(3.17) 

We write 


Vi (zo) + V(y, x 0 ) V1 + d(y, x 0 ) 

1 


X{y£X-.d(y,x 0 )>R/8}'-i 


'X 


\h{y)\dn{y) < CWfUM 


1 + d(x, Xq) J 


h{v) dfx(y) 


'X 


fi(y)dy(y)~ / f 3 (y) dfi(y) 


ix 


< CWfWwrt 


J x 

1 \7 


1 + d(x,x Q ) 

T(fi)(x) = I K(x,y)u(y)[fi(y) — fi(x)]d{j,(y) 


'X 


P(x) 

K(x,y)v(y)f 1 (y)d^(y) +f^x) / K(x,y)u(y) d[i(y), 


lx 


lx 


-V- 

'IX! 


where u(y) is a smooth cut-off function as in Lemma [3781 satisfying u(y) = 1 when y e B(x,r), 
and u(y) = 0 when y G B(x, 4Agr) c , and where v(y) := 1 — u(y). 
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Iii order to estimate the expressions p(x) and q(x), we show that the kernel K(x,y ) of T 
satisfies the following four estimates: 

C 


(3.18) 

for all x f- y\ 

(3.19) 

for d(x,x t ) < ( 2 A 0 ) _ 1 d(:r, y)‘ 


\K(x,y)\ < 


V(x,y) 


I K(x,y) - K{x\y)\ < ° ^ j Y 

V{x,y) V d{x,y) J 


(3.20) 


\K(x,y)-K(x,y')\ < 


C 


( d(y,y') V 


V(x,y)\d(x,y) ) 

for d(y,y') < (2 A 0 )~ 1 d(x,y)] and 

(3.21) |*(*,») - K ix ,y>) - K M ) + K {x ',y ')I < 

for d(x,x') < ( 2 A 0 ) _ 1 <i(;r, y) and d(y,y') < {2A 0 )~ l d(x,y). 

We point out that without additional assumptions on the measure /i, the decay and smooth¬ 
ness estimates (l3.9jMI3.10ll for D k (x,y) as given in Lemma [3T6l are not by themselves suf¬ 
ficient to imply the above estimates for K(x,y ) = Dk(x, y). Fortunately, however, in 
our setting, although we have no additional assumptions on /./, we do have the special form 
K(x,y) = ZkDt (x,y) = Y7k a ' l l J a( x )' l l J a(y) i n terms of the wavelet basis, rather than the 
operators D k = S k +1 — S k from the classical case. Instead of using the estimates (13.911 - 
(13. 10H directly, we can use the estimate (I3.22jl proved below, together with the approach of 
Lemma 9.1 in EH]. The estimates (13.18jl — (13.2011 for our K(x,y) were proved in Lemmas 9.2 
and 9.3 of (sg. Thus, we only need to show the estimate (13.2111 . To do this, following 
the approach of Lemma 9.3 in eh, we first claim that if d(x,x') < (2A 0 ) 1 d(x,y) and 
d(y,y') < {2A 0 )- l d(x 1 y) 1 then 


(3.22) 


E | KW - - €{y')\ 

C . r. / , d{x,x')\ ri 


< 


V(x, 5 k ) 


mm 


5 k 


nun 


d(y,y ')v 
S k 


x exp(— k d(x, < 3f k )) a ) exp(— v{8 k d(x,y)) a ). 

Recall that a := (1 + log 2 A 0 ) _1 is the exponent defined in (j2.9j) in Theorem 12.21 

We prove (13.2211 following the method used to prove the second assertion in Lemma 9.1 
of [AH]. We consider four cases. First suppose 5 k > d(x,x') and 5 k > d(y,y'). Then 

[V£(®) - ^a( X ')\ bPaiv) - ^a(y ')] 

'd(x,x')\v fd(y,y')\v 


< 


< 


C 

V(x, 5 k ) 
C 


5 k 


5 k 


V (x, S k ) \ S k 


d(x,x')\ J }/d(y,y')\ r > 


5 k 


eM-v(d- k d(x, y k )) a ) exp (-v(8~ k d(y, y k )) a ) 


exp (~u(5 d(x,y a )) a )exp(-v(5 ~d(x,y)) a ). 


(Here the constant v changes from line to line to accommodate the constant Aq that arises 
from the use of the quasi-triangle inequality.) The sum over a G *3/ k of the fourth factor on 
the right-hand side is dominated by the expression exp(— u(5~ k d( y x 1 3^ k )) a ). 
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Second, suppose 5 k > d(x,x') and 5 k < d(y,y') < (2 A 0 ) 1 d(x,y). Then we have 


M(*)-V£(*0] K(v) ~ €(y')] 


< 


bl>a(x) - t/>aW)] Ipaiv) + “ ^a(^)] V£(Z/') 


Then, by using the second estimate in Lemma 9.1 from [AH] , and the quasi-triangle inequality, 
we obtain that 


KW - [€(v) - €W)} 

< (-^T —yexp(-v(5~ k d(x, yj))“) exp (-v(5~ k d(x t y)) a ). 

Here the sum over a G ( 3/ k of the third factor on the right-hand side is dominated by 
exp (—u(8~ k d(x, k )) a ). 

The other two cases, namely when 5 k < d(x,x') < (2A 0 )~ 1 d(x, y) and 5 k > d(y,y'), and 
when 5 k < d(x,x') < (2 A 0 )^ 1 d(x,y) and S k < d(y,y') < (2A 0 )~ 1 d(x,y), can be handled 
similarly. We omit the details. 

Combining the estimates for all four cases above, we have established the claim (13.221) . 
Now we verify (13.211) . From the definition of K(x,y ) and the claim (13.221) . we see that 


I K(x, y) - K(x, y ') - K(x, y') + K(x’, y’) \ 

< | WO*) _ ^a( x ')\ bP*(y) - V£(s/')] 

k a£^ k 

S E 

k: 5 k >(2Ao) 1 d(x,y) 

+ E 

k: 5 k <(2A 0 )~ 1 d(x,y) 


+ E 

k : S k <(2Ao) 1 d(x,y) 
d(x,x')>8 k ,d(y,y')<5 h 

+ E 

k: S k <(2A 0 ) _ 1 d(x,y) 


+ 


E 


k: 8 k <(2Ao) 1 d(x,y) 
d(x,x')>S k ,d(y,y')>S k 

—■ B\ + B 2 + B 3 + B4 + H5. 


c ( 

<d(x, x') 

Y(d(y,y')\ 

V 

1 exp(— i/(S' 

V(x, S k )' 

v 5 k 

) V a* J 

C | 

(d(x,y)' 

fd(x, x f )\ 

v (d{y,y')\ 

v(x, y y 

< 5 k / 

1 \ 8 k ) 

V 5 k ) 

c | 

( d(x, y) ' 

\“(d(y,y')\ 

V 

exp(— u(5" 

V(x,y )' 

^ S k y 

I \ 5 k ) 

C | 

( d(x,y)' 

fd(x, x')\ 

V 

exp(— u(S' 

v{x, y y 

< 8 k y 

1 \ 8 k ) 

c | 

(d(x,y)' 

1 exp(— u(5 

~ k d(x, y)) a ) 

v{x, y y 

< 8 k y 


k d(x, & k )) a ) 
exp(-z ;(5- k d(x,y)) a ) 


Following Lemma 8.3 in EE], an application of the estimate in Remark 12.31 shows that 
B\ is bounded by v ^. y ^ ( d j^ y j ) V ( ‘d^y^ Y’ ■ Further, B 2 satisfies the same estimate since 
Ylm= 0 5 _m ^ +T? ) exp {—u5~ ma } < C. We can deal with B 3 , R 4 and B 5 similarly. Thus (13.211) 
is proved. 

We remark that the estimate (I3.21[) is crucial for the proof of (13.12ji : see the estimate for 
T(f 2 ) below. 
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Now that we have established the estimates fl3.18jl - fl3.21h on the kernel K(x,y) of T, we 
return to estimating the expressions p(x) and q(x). The size condition on the kernel K(x,y ) 
and the smoothness condition (13.1411 on /i yield 

IpWI < cimi 

1 \ 7 

dy(y) 


x 


hi(x 0 ) + V(x, x Q ) V1 + d(x, x 0 ) 


/ d(x,x') \ 

f 1 ( 1 ] 

Vl + d(x, xq) J 

Vi (xo) + V(x, x 0 ) V 1 + d(x, x 0 ) ) 


This estimate still holds when x is replaced by x', for d(x, x') = r. Thus 

| PW - P MI < cil/llo, < v,)( 1 + V( X , X0 ) ( i + dL,) ) 7 - 

For q(x), since T1 — 0 (by the definition of Dk(x,y) and the cancellation property of ip*), 
we obtain 

q(x)-q{x') = [ [K(x,y) - K(x',y)]v(y) [f^y) -/i(x)] d/x(j/) 


I x 


+ [h{y) ~ h(x)] / K(x,y)u{y)dy(y) 
=: (E) + (F). 

We claim that there exists a constant C such that for all x, 


(3.23) 


K(x, y)u(y) dy(y) 


i x 


< C. 


Assuming this claim (which is proved below), together with the estimate for fi in (13.1411 . we 
hnd that 

|(F)| < cimx) -/.Ml < CI/l| c „,,(A||T y ) % i(io) + V feio) (^- 

Applying the smoothness estimates for both f\ and K(x,y), we obtain 


1(E) | < C 


d(x,y)>4AQr 


I K{x,y) - K(x',y)\ |u(r/)| |/i(?/) - /i(x)| dy(y) 

1 


< CII/lloOS.T) V(i, y) ( d(x, y) ) 


X 


( d(x,y ) \ 

{ 1 | 

( 1 \ 

Vl 4- d(x,xo)J 

Vi(x 0 ) + V(x,x 0 ) 

Vl + d(x,x o)J 


d(x,x') 


dy(y) 


Y 1 i 

< 1 1 

' Ui(x 0 ) + U(x,x 0 ) 

Vl + d(x,x 0 )J 

d(x,x') 

1 ( 

l + d(x,x 0 )J Vi(x 0 ) + V(x, x 0 ) Vl 


since f3 < rj. Therefore 


We consider three cases. First suppose d(x,x') = r < (20Aq) *(1 + R) and R > 10. Then 
the points x and x' are not in the supports of / 2 and f?,. Using the double smoothness and 
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smoothness conditions ( (13.211) and (13.19j) respectively) on K(x,y ), and the estimate (13.171) 
of / 2 , we find 


\T{f 2 ){x)-T{f 2 ){x')\ = 


\K(x,y ) - K(x',y)]f 2 (y) dfi(y) 


' x 


< / \K(x,y)-K(a?,y)-K(x,xo)-K(rf,xo)\\f 2 (y)\d[i(y) 


>x 


+ \K(x,x 0 ) - K(x',x 0 )| 

<cmwJ d{x ’ x ' 


i x 


h{y) dy(y) 
1 


1 + d(x,x 0 )J Vi(x 0 ) + V(x,x 0 ) Vl + d(x,x 0 ) 


Also, 


\T(f 3 )(x) ~T(f 3 )(x')\ = 


[K(x, y) - K(x', y)\h(y) dfi(y) 


ix 


< C 


' d(x ,y)> Q>2Aor 




<C\\f\\QU,„) 


d(x, x') 


.1 + d(x,x 0 )J Vi(x 0 ) + V{x,x Q ) Vl + d(x,x 0 ). 

In the second case, where d(x,x 0 ) = R and (2A 0 ) _1 (l + i?) > d(x,x') — r > (20Aq) _1 (1 + .R), 
the desired estimate for T(f)(x) follows from the estimate of (13.61) . So we need only consider 
the third case, where R < 10 and r < 11/(20A 0 ). This case is similar, and indeed easier. In 
fact, all we need to do is to replace R in the proof above by 10. We leave the details to the 
reader. This completes the proof of (13.71) . 

To finish the argument for Theorem 13.41 it remains to establish the claim (13.231) . To do 
so, we prove that there exists a constant C such that 


(3.24) HT0IU < C 

for all functions 0 with the properties that ||0||oo < 1 and there exist x 0 G A" and t > 0 such 
that supp 0 C B(x 0 ,t ) and ||0||, ? := sup x ^ y {\(p(x) — (f>(y)\/d(x : y) v } < t~ v . 

We again follow the idea of Meyer’s proof in |M2j . Let Xo( a; ) — h(x), where h(x) is a 
smooth cut-off function as in Lemma [3.81 with the property that h(x) = 1 on B(xo,2t) and 
h(x) = 0 on B(x o, 8Alt) c . Set xi := 1 — Xo- Then 0 = 0y o and for all 0 G Cq(X), 

(^ 0 , 0 ) = (K(x,y),<f>(y)il)(x)) = (K(x,y),xo(y)(/>(y)^(,x)) 

= (K(x,y),Xo{y)[<f>(y) - 0(aO]0(z)) + (K(x,y),Xo(y)<f>(x)iJ;(x)) 

:= (G) + (H). 

Applying the size condition (13.181) on the kernel K(x,y ) yields 

l(G)|<C||0|| 1 . 

To estimate (H), it suffices to show that for x G B(x 0 ,t), 

(3.25) \T X o(x)\<C, 
since as (H) = (Ty o ,00), we then have 

I (H) | < 11 Ty o 11 ijco (b(xo ,*)) 11 faP IIl 1 (b(xq ,t)) <G||0||i. 
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To show (13.251) . we use Meyer’s idea again [M2]. Take ip G C\X) with supp-0 Q B(xo,t) 
and f x ip(x) d/i(x) = 0. Since T 1 = 0 and f x ip(x) dfj,(x) = 0, and using the smoothness 
condition (13. 19ft on K(x,y), we obtain 


\{ T Xo,i>)\ = \~( T Xi,^)\ = 

< emu. 


[K(x, y ) - K(x o, y)]xi(y)^(x) dfi(x) dfi(y) 


XxX 


Thus, T\o(x) = A + j(x) for x G B(x 0 ,t), where A is a constant and ||7||oo < C. To 
estimate A, choose (pi G Cq(X) with supp^i C B(xo,t), ||0i||oo < 1, ||0i||»? < t~ v and 
f x 4>i(x) dfi{x) = Ct. Since T is bounded on L 2 (X), we have 


CtA + 


(pi(x) n f(x) dy(x) 


\(Txo,<Pi) \ < Ct. 


J x 

Therefore |A| < C, and hence the claim (13.23ft is proved. This completes the proof of 
Theorem 13.41 modulo the proof of Lemma 13.61 □ 


It remains to prove the technical lemma used in the proof of Theorem 13.41 


Proof of Lemma Id. hi (i) To establish the decay condition 


, we write 
^a(y) 






By Theorem 13.31 we know that 'ip k ( x )/\Jy,(B(y k , S k )) belongs to G(y k , 6 k , rj, 7 + (]). Applying 
the size condition (i) from Definition 13.11 we see that 

\o t (x,y)\ < C £ M (s(rf,^)) VM) + v(yJix) + dbA ' X) . 




X 


7+r? 


V S k{y*) + V(y k ,y) + d{y k ,y). 

Note that for each 2 G B(y k , d k ) one can replace y k by z to get 5 k + d{y k , x) ~ 5 k + d(z, x) 
and V S k(y k ) + V(y k ,x) ~ V S k{z)+V{z, x), and similarly for S k +d(y k , y) and V S k(y k )+V(y k , y). 
Thus, first replacing y(B(y k ,S k )) by j R(y k S k) dy(z) and then replacing y k by z, and finally 
summing up over a G l 3 / ' k , we find that the last sum above is bounded by 


C 


( 1 

/ 5 k 

\ z+v 1 / 5 k \ 

Jx V S k (z) + V (z, x) 

\5 k + d(z,x). 

/ V S k(z) + V(z,y)\S k + d(z,y)J 


7+7? 


d/j,(z) 


=■ (P) + (Q). 


Here (P) is the result of integrating over the set d(x, z) < ( 2A 0 )~ 1 (5 k + d(x, y)) and (Q) over 
the set d(x,z) > (2A 0 ) _ 1 (h fc + d(x,y)). To estimate (P), note that if d(x,z) < (2A 0 ) _ 1 (h fc + 
d(x, y)) and 2 5 k < d(x, y), then d(y, z) > ( 10 A 0 )^ 1 (h fc + d(a;, y)) and by the doubling property, 

V(z,y) = fi(B(y,d(z,y))) > y( B (y, ( 10 A 0 ) _ 1 d(x, y))) > ( 10 A 0 )-“V(x,y). 


Therefore, V S k(x) + V(x,y) < C(V 5 k(z) + V(z,y)). Next, if d(x,z) < (2A 0 )~ 1 (S k + d(x,y)) and 
2 5 k > d(x,y), then d(x,z) < 3(2 A 0 )~ 1 5 k . For this case, first suppose that d(x,z) < 5 k and 
hence, V S k(z) ~ V S k{x). On the other hand, if d(x,z ) > 5 fc , then (2A 0 —l)<5 fc < (2 A 0 )^ 1 d(x,y) 
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and hence V 5 k(x) < CV(x,y). Therefore, in this case, again V S k(x ) + V(x,y) < C(V S k(z) + 
V(z,y)) and thus we get 


(P)<C 


V S k (x) + V(x, y)\S k + d(x, y) 


5 k 


7+4 


< C 


V S k (x) + V (x, y) V 5 k + d(x, y) 


S k 


7 

? 


as required. The estimate for (Q) is the same, but with x and y reversed, 
(ii) To establish the smoothness condition (13.101) . we write 


Dk(x,y) - D(x,y')\ 

= E 


^0 


X 


^a(y) 


V^( B (y*’ 5k )) W^ B (ya,^ k )) 


=■ (R) + (S). 


^a(y') 

Vv( B (ya>ti k ))- 


ffere (R) is the result of summing over the set of a G , 3^ k such that d(y,y') < (2Aq) 1 (S k + 
d(y,y k )) or d(y,y ') < (2y4 0 ) _1 (^ fc + d(y',y k )), and (S) over the set of a 6 r 3d k such that 
d(y,y') > (2A 0 )-\S k + d(y,y k )) and d(y,y') > (2 A 0 )~\5 k + d(y\ y k )). 

For (R), use the size condition (Definition 13.If i)) for the first factor ift k (x )/ y/ y(B(y k , 5 k )) 
and the Holder regularity condition (Definition 13.1( h)) for the terms ^ k (y)/\fy,(B(y k , 6 k )) 
in the second factor. We find that 


(«)<E y(B(y k a ,6 k )) 

<X&8f k 

1 


( d(y, y') \ 

v 1 / 5 k \ 

\5 k + d(y k , y)) 

1 v 5 k (y k ) + V(y k ,y) + d(y k ,y)J 


7+4 


X 


Vgk (y k ) + V (: y k , x) \S k + d(y k , x) 


\ 7+4 


Applying the same proof as for (13.91) . we see that the last sum above is bounded by 


f f d(y, y') y 

1 , 

f S k \ 7+4 

lx \S k + d(z,y)J 

v 5 k(z) + v{ z , y y 

^8 k + d(z, y) ) 


1 ( d k y+v 

V$k(z) + V(z,x)\5 k + d(z,x)J 


< r( d ( y '^ V 1 ( 6k V 

— \5 k + d(x, y)J V S k (x) + V(x,y)\5 k + d(x,y)/ 


To deal with (S), we can write 


(S) < E k-{B(vi,S k )) 

ae@ rk :d(y,y')>(2Ao)- 1 (6 k +d(y,y k )) 


^ k (x) 


^a(y) 


+ Y, ti B (Va,8 k )) 

a&'3f k -.d(y,y')>(2Ao)- 1 (5 k +d(y l ,y k )) 


^fi(B(y k ,5 k )) a/ d( B (y k , 5 k )) 

V^O) iPa(y') 




For the first sum, following the same approach as for (R) but with d(y,y') > (2 A 0 ) 1 (S k + 
d(y,y k )), we must deal with the integral 


f ( d(y, y') y 

1 | 

( s k \ 7+4 

1 | 

( s k \ 

J x \5 k + d(z, y)) 

V S k(z) + V(z,yy 

^ 5 k + d(z, y)J 

V S k(z) + V(z,x) 

^5 k + d{z,x)J 


7+4 


dy(z). 
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Applying the same proof as for (R), but using the size condition (Definition 13.If i)) for both 
factors, we obtain that this integral is bounded by 

c ( d(y,y') y 1 / s k y 

\S k + d(x,y)J V S k(x) + V(x,y) \8 k + d(x,y)) 

The second sum is similar to the first one, with y and y' reversed. Thus, by the same proof 
we find that the second sum is bounded by 


( d(y, y') \ 

11 1 

( 6k \ 

\8 k + d(x,y')) 

V 5 k(x) + V{x,y') 1 

^5 k + d(x, y') J 


Note that the fact d(y,y') < (2A 0 ) _1 (5 fe + d(x,y)) implies that 8 k + d(x,y) ~ 5 k + d(x,y') 
and V 5 k(x) + V(x,y) ~ V S k(x) + V(x,y'). Therefore, we obtain the desired estimate for the 
second sum. 

(iii) The proof for the double smoothness condition (13. lip is similar to that for (13.101) . and 
we omit the details. 

This completes the proof of Lemma 13.61 □ 

3.2. Product test functions, distributions, and wavelet reproducing formula. We 

now consider the product setting (Xi,di,pbi) x (X 2 , d 2 , 7^2), where (A 7 ,, di, 77), i — 1, 2, are 
spaces of homogeneous type as defined in the Introduction. For i — 1, 2, let be the 
doubling constant as in inequality (jl . 2[) . let 07 be the upper dimension as in inequality (11.31) . 
and let be the constant in the quasi-triangle inequality (II.ip . In this subsection we use 
the notation (x,y) for an element of X\ x A 2 . 

On each X\ there is a wavelet basis with Holder exponent ry as in inequality (12. lip . 

We now define the spaces of test functions and distributions on the product space Ah x A 2 . 


Definition 3.9. (Product test functions) Let (xo,yo) £ Ah x X 2 and r = (rr,r 2 ) with ri, 
r 2 > 0. Take /3 = (/h.,/h), with /h G (0,771], /3 2 £ (0,772], and 7 = (71,72) with 71, 72 > 0. A 
function f(x, y) defined on Ah x A" 2 is said to be a test function of type (xo, t/o) D /3] 7) if the 
following three conditions hold. 


(a) For each fixed y £ X 2 , f(x,y) as a function of the variable x £ Ah is a test function 
in G(ao,ri,/7i,7i). 

(b) For each fixed x £ Ah, f(x,y ) as a function of the variable y G A" 2 is a test function 
in G(yo, r 2 , /3 2 ,72). 

(c) The following properties hold: 

(i) (Size condition) For all y G X 2 , 

1 f r 2 \ 72 

"/(-.!/)Ikw..-n> < C VM + v(v m) ( r2 + ■ 

(ii) (Holder regularity condition) For all y, y' G A 2 such that d 2 (y, y') < (2AQ 2 ^ 1 (r 2 + 
d 2 (y,yo)), we have 




xo,n,/ 8 i, 7 i) — 


< c 


d2(y,y') 


r 2 


72 


-r 2 + d 2 (y,yo)J V r2 (y 0 ) + V(y 0 ,y) \r 2 + d 2 (y,y 0 )- 

(iii) Properties (i) and (ii) also hold with x and y interchanged. 

(iv) (Cancellation condition) f x f(x,y) dy,i(x) = 0 for all y G Ah, and / V2 f(x, y) 
dy 2 (y) = 0 for all x G Ah- 
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When / is a test function of type (#o, Vo', r ", ft 7), we write / G Gfto, 7/0; r; ft 7). Note the 
use of semicolons here to distinguish the product definition from the one-parameter version. 
The expression 

ll/llG(*o,w;r;/»;7) : = ■ (i), (h) and (iii) hold} 

defines a norm on G(x 0 , y 0 ] r; ft 7). 

We denote by G(/ 3 1, ft; 71,72) the class G(x 0 , y 0 ; 1 , 1 ; ft 7) for arbitrary fixed (x 0 ,y 0 ) G 
x X 2 . Then G(x 0 , y 0 ] r; ft 7) = G(ft, ft; 71,72), with equivalent norms, for all (x 0 ,y 0 ) G 
Ab x X 2 and ri, r 2 > 0 . Furthermore, G(ft, ft; 71,72) is a Banach space with respect to the 
norm on G(/ 3 1,/ 3 2 ;71, 7 2 )- 

o 

For Pi G (0,7/j] and 7* > 0 , for i — 1 , 2 , let G(ft,ft;7i,72) be the completion of the space 
G(r) 1,72; 71,72) hi G(ft,ft; 71,72) in the norm of G(ft, ft; 7i, 7 2 ). We dehne the norm on 

G(ft, ft; 7 i) 72) by II/II S(/JiA;71:72) := ll/llG(A,fe ;7l , 72 )- 

The (scaled) product wavelets given by 1/2 

are product test functions in G(y^\, y^ 2 ; 8 ] ft 7) for each 7 = (71,72) with 71 > 0, 72 > 0, 
where 5 = (d}’ 1 , d/? 2 ) and ft = [r] 1,772); this is straightforward to check. 

Definition 3.10. (Product distributions) Let (x 0 ,yo) G X 1 x A " 2 and r = (r!,r 2 ) with 77 , 

r 2 > 0 . Take ft = (ft, ft), with ft G (0,71], ft G (0,7/2], and 7 = (71,72) with 7!, 72 > 0 . 

We define the distribution space (G(Pi, ft; 71,72)/ to consist of all linear functionals C from 
0 _ 

G(ft, ft;71,72) to C with the property that there exists a constant C such that for all 
/ G G(ft,ft; 71,72), 

We have the following version of the wavelet reproducing formula in the product setting 
Ab x X 2 . 


Theorem 3.11. (Product reproducing formula) Take Pi, 7 \ G (0,7ft /or* = 1 , 2 . 

(a) The wavelet reproducing formula 

( 3 . 26 ) f( x , y ) = J 2 

fcl fc 2 ot 2 &<3f k 2 

O 

holds in the space of test functions G(ft,ft; 7i,7 2 ) / or eac/i Pi G ( 0 , ft) and 7' G 
( 0 , 7i), fori = 1 , 2 . 

(b) Tde wavelet reproducing formula (I 3 . 26 jl also holds in the space of distributions 

(G(ft,ft; 71,72))'- 

Proof. As before, the wavelet reproducing formula for distributions follows immediately from 
that for test functions. The proof for test functions proceeds by iteration of Theorem 13.41 
Write 

9 (x,y) := E E E E f( X ^y) 

\ki\<L\ ai£:3f k i \k 2 \<L 2 a 2 £?V k 2 

=■■ 9 ii x ,y) + 92(2,9), 
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where 

9 i(x,y) 


and 


E E (W E E (/(■•■),}<w 

\ki\<Li \k 2 \<L 2 a2 6^ fc 2 

“EE (/Ow), ^aD^aliv) 

\& 2 \<L 2 Q2€^ fc 2 


92 (x,y):= E E (/( x u),^S)^a 2 (2/) ~ f(x,y). 

\&2\<L2 «2G^ fc 2 


To see the convergence in the space of test functions, we recall the following (one-parameter) 
estimate on A", as shown in the proof of Theorem 13.41 Given /?, 7 G ( 0 , 7), for each f 3 ' G ( 0 , ( 3 ) 
and 7' G (0,7) there is a constant a > 0 such that for each positive integer L 


( 3 . 27 ) 


/(•)- E Eftwso 


|fc|<I/aG^ fc 




GCS',7') 


where C is a constant independent of / G G(/ 3 , 7). Note that inequality (I 3 . 27 jl is the same as 
inequality ( 13 . 5 ) 1 . slightly rewritten. Inequality ( 13 . 27 )) . together with the triangle inequality, 
implies that 


( 3 . 28 ) 


E EW’.v)^(•) 

|fc|<L 


G(/3', 7 ') 


< C\\f\\ov„). 


We observe that if / G G(f 3 1, / 3 2 ;7i,72), then ||/(-,2 /)||g(,Si,7i), as a function of the vari¬ 
able y, is in G(/ 3 2 , 72), and satisfies ||||/(-, •)IIg(/ 3 i, 7 i)|| g(/ 3 2)72) - ll/ll<309i,/3a;7i.72)■ Similarly, 
|| !!/(•> OIIgJ/Sh, 72) IIg^.tO ^ II/IIg(/3i,/3 2 ;7i,72)- Therefore, we obtain 


lG(/3i,7i) 


||5'i(-,2/)||g(^, 70 < ^ Ll 1l E E </(■»■),^S(y)|| ( 

1^21< 7^2 a 2 e^ fc 2 

■ )l| o^)i^rn^,y) Vr 2 + %,?/o)^ " G{M 


< CA Lic 


r 2 


72 . 


< CV> Ll °'||/(-,-) ||G(/ 3 i,/ 3 2 171,72): 


P 2 


\ 72 




K 2 (z/o) + v(y 0 ,y) Vr 2 + d(y,y 0 ). 

where the first inequality follows from ( 13 . 271 ) and the second inequality follows from ( 13 . 281 ) . 
Similarly, 


Il02(a:, Z/)llG(/3i,7i) < C5 L2a 11 / 11 G(fh ,/9 2 ;7i > 72 ) 


Ka(z/o) + V(yo,y) Vr 2 + d(y,y 0 ) 


^2 


72 


Noting that g(x,y) - g{x,y') = \gi(x,y) - gi(x,y')] + [g 2 (x,y) - g 2 (x,y')\, by repeating the 
same estimates we obtain 


\\g(-,y) - gi^y^Wawri) 

<C(S L ^ + 5 L n\\f(;-) 110(^2171,72) 


( d(y, y') \ 

.* 1 1 

f ) 

\r 2 + d(y,y 0 )) 

Vr 2 (yo) + V(y 0 ,y) 

W 2 + d(y, yo)) 


where d{y,y') < (2A ( 0 2) ) 1 (r 2 + d(y,y 0 )). 


72 
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The same proof can be carried out for the estimates with x and y interchanged. Hence 
IM|g(/3P/3'; 7 p 7 ') < C(S L ^ + S L n\\f\\ G{h ,y m h 

o 

which yields the convergence in G(P[, fy] 7i, 72 )- □ 

4. Littlewood-Paley square functions and Plancherel-Polya inequalities 

We now carry out the philosophy described near the end of the introduction, in order to 
establish the Littlewood-Paley theory for the discrete square function in terms of wavelet 
coefficients. We define the discrete and continuous square functions, and prove their norm- 
equivalence via Plancherel-Polya inequalities, whose proof takes up most of this section. 
Again we begin with the one-parameter case. 

4.1. One-parameter square functions via wavelets, and Plancherel—Polya inequal¬ 
ities. We first apply the orthonormal wavelet basis constructed in [ AH ] to introduce the 
discrete Littlewood-Paley square function, defined via the wavelet coefficients as follows. 

o 

Definition 4.1. (Discrete square function in terms of wavelet coefficients) For / in (G(/3, 7 ))' 
with / 3 , 7 € ( 0 , 7 ), the discrete Littlewood-Paley square function S(f ) of f is defined by 

(4- 1 ) S(f){x) := KV&/>XQfc(aO| 2 } ‘ , 

k a £<y k 

where \Qk (x) := XQ k ( x )h(Qa)~^ 2 an d XQ k ( x ) is the indicator function of the dyadic cube Q k a . 

It is straightforward that \\S(f)\\ L 2 ^ X ) = || /1| z. 2 (at) 5 since { 1 } forms an orthonormal wavelet 
basis for L 2 (X). However, it is not easy to see why ||<S , (/)||lj>(.x;) rs j \\f\\ L p(x) for 1 < p < 
00 with p ^ 2. This difficulty is because the classical method, namely the vector-valued 
Calderon-Zygmund operator theory, cannot be carried out here due to the lack of smoothness 
in the x variable. For this reason, we introduce the following continuous Littlewood-Paley 
square function in terms of the wavelet operators Dk- 

Definition 4.2. (Continuous square function in terms of wavelet operators) Let Dk be the 

O 

operator with kernel D k (x,y) = Y^ a e& k ^a{ x )$a(y)- For / e (G{/3, 7))' with ^,76 ( 0 , 7 ), 
the continuous Littlewood-Paley square function S c (f ) of f is defined by 

Sc(f)(x) := \D k (f)(x )\ 2 } ' . 

k 

The two main results in this subsection are as follows. 

Theorem 4.3. (Littlewood-Paley theory) Suppose (3, 7 G ( 0 , 7 ) and < p < 00 , where u 

o 

is the upper dimension of (X,d, p). For f in (G(/3, 7))', we have 

\W)\\ LP(X) ~ \\s c (f)\\ Lnx) . 

Moreover, if 1 < p < 00 , then 

IW)|UqX) ~ ||5 'c(/)||lp(X) ~ \\f\\LP(X)- 

The key idea in proving Theorem 14.31 is the following Plancherel-Polya type inequalities. 
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Theorem 4.4. (Plancherel-Polya inequalities) Suppose ^76 (0, 77 ) and < p < 00 , 
where u> is the upper dimension of (A", d, p). Fix N G N. Then there is a positive constant C 

o 

such that for all f G (G(/3, 7 ))', we have 


{E E [ SU P \ D k'{f)(z)\ 2 ]x Q k'+N(-)\ 

1 " l ZZQ k f+ N J 

E E I (^f)xQ^) 


1/2 


Lp(X) 


< C 


1/2 


Lp(X) 


(42) 

k a e9' k 

Moreover, for a fixed sufficiently large integer N (N will be determined later in the proof), 

o 

there is a positive constant C such that for all f G (G(P, 7 )/, we have 


E E 

k o.e?y k 


1/2 


LP(X) 


(4.3) 


< C 


k ' a'&5£ k '+ N 



,1 1/2 

X q y+ N {- 

J Ct' 



LP(X) 


Note that in each of the inequalities (14,2ft and (14.31) . on one side, for each k G Z the sum 
runs over the set & k , while on the other side for each k! G Z the sum runs over the set 3F k '+ N . 
Besides the distinction between and SF, the other difference here is that in the expressions 
involving D k *, it is not sufficient to sum at the scale of k', but rather, following ra , we 
must sum over all cubes at the smaller scale k' + N. 


Proof of Theorem \4--i\ Theorem 14.31 follows from Theorem 14.41 by standard arguments that 
can be found in EH]. We sketch the idea. The first estimate in Theorem 14.31 follows from 
Theorem 14.41 together with the following observation: 




E 


inf 

k' a 'ex k '+ N zeQ * 


k' + N 


\D k/ (f)(z)\ 2 x Q k'+ N (x) < El Dk(f)(x)\* 


< 


E 


x sup \Dk>(f)(z)\ x k' +N fx). 

a'e5l~ k '+ N 2 £< 9 a' +JV 




For the second estimate in Theorem 14.31 when 1 < p < 00 one obtains from the classical 
method of vector-valued Calderon-Zygmund operator theory that ||5 , c (/)||z /P (x) < C||/||lp(x)- 
This estimate together with the wavelet expansion as in (12.13j) gives ||/||i>(x) < C\\S c (f)\\LP(x), 
and Theorem 14.31 follows. □ 


We would like to point out that to consider S c (f) as a vector-valued Calderon-Zygmund 
operator, we need to use the crucial estimate mentioned in Remark l2.3l to show that the kernel 
of the operator S c (f) satisfies all conditions for the Calderon-Zygmund singular integral 
operator. We omit the details. 

Outline of proof of Theorem 14.41 Since the proof (below) is rather complex, we begin 
by outlining our approach. For the first Plancherel-Polya. inequality (14.21) . we substitute the 
wavelet reproducing formula (13.31) for / into the left-hand side. Thus the desired wavelet 
coefficients (if^, f) appear. To deal with the unwanted terms D k ' and we apply the 
almost-orthogonality estimates (14.41) given below. Then the standard technique, as in 1)11,. 
of applying an estimate from Ea and the Fefferman-Stein vector-valued maximal function 
inequality [ FS] establishes (14.21 . 
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The second Plancherel-Polya inequality (14.3j) is harder. Roughly speaking, we need to 
control the wavelet coefficients by the quantities D k '{f). Now for spaces of homogeneous 
type with additional assumptions, one proceeds as in [}DH] via a frame reproducing formula 
of the form 

/(*) = E E ^Q^ +N )D k '(x,x k ^ +N )D k (f)(x^ +N ). 

k' a'&3C k '+ N 

However, for our spaces of homogeneous type with no additional assumptions on d and //, no 
such frame reproducing formula is available. A new idea is needed. We introduce a suitable 
operator Tjy, show that X n is bounded and that the L P (X) norm of S , (T / ^ 1 (/)) is controlled 
by that of S(f ) (Lemma 14.61 below), and rewrite the wavelet coefficient as ('ijj k ,T^ 1 T N f). 
Pulling out the operator T^ 1 from the left-hand side of (14.31) . we obtain expressions of 
the form {ip*/\Jp(Qa)i T N f). Because of the form of T N , we can now apply the almost- 
orthogonality estimates (14. 4 j) to these terms and complete the remainder of the proof of the 
second Plancherel-Polya inequality (14.31) by following the approach used for (14. 2 j) . 

We now give the details. 


Proof of Theorem \4. f\ For each / G {G{/3, q))', by Theorem 13.41 the functions 

•Mi 


fn(x) = E E 

\k\<na&W k 


\X) 


belong to L 2 {X) and converge to / in (G(/3, 7 ))' as n —> 00 . Note that (/„, ip*) = {f,ip k ) for 
| A:| < n, and (f n , ip*) = 0 for \k\ > n. Thus, 

12 " 


E E sup \D k ff)(z)\ x Q k'+N(x) 


\k'\<n a'&3(' k '+ N z€ Qt ,+N 


E E SU P \ Dk '(fn)(z )\ 2 Xn k ’+N{x) 


k’ a ' &3( :k'+N " z£Q k f +N 


and 


E E K ^f)XQ k 

|/c|<n 


E E I (^aJn)XQ k 

k aeffc 


Therefore it suffices to show the inequality (14.21) of Theorem 14.41 for / G L 2 {X), and similarly 
for the inequality (j4.3[) . 

We first prove (14.21) . Fix iV G N. The idea is to apply an almost-orthogonality estimate 
( (14.41) below). First, for each / G L 2 {X), by the wavelet expansion (Theorem 13.41) . 

/w = X Ed’Ww- 

fceZ aG&' k 


Thus for each z G Q k : , +N we have 


\/ ^a(') 




Claim: ( Almost-orthogonality estimate ) We claim that ('ip k (')/VP(Qa)^ Ac'G, z)) satisfies 
the following almost-orthogonality estimate: There exists a constant C such that for each 
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positive integer N, each 7 > 0, each point z G Q^' +N and each point x k a A N G Q^ +N , we have 

(4.4) < 


§(kAk') 


»W*S) + V*>(* 5 +,r ) + v(*t>*Z+ H ) Vi5( ' sM ' ) + 

As usual, k A k 1 = minjfc, A;'} denotes the minimum of k and A/. 

Remark 4.5. The key idea used below to prove the claim (j4.4[) is that both i^a( x )/V^(Qa) 
and D k /(-,z) satisfy size conditions, Holder regularity conditions, and cancellation, since 

as we have shown, '0«( a; )/'\// i ( ( 3a) * s a t es t f ,inc tion in G{/3, 7 ) while D k >(-,z) satishes the 
properties fj3.9jl — (13. 1 ljl in Lemma 13.61 Further, we point out that if D k (x,y ) satishes the 
same size condition (13.9|) together with the following Holder regularity condition (which is 
weaker than (13.1 OK ). 

'd(y,y')\v r 1 / S k \7 


\D k (x,y) - D k (x,y')\ < C 


5 k 


L V S k (x) + V(x, y)\S k + d(x, y) 


5 k 


71 


^ ^ V 5 k(x) + V(x,y') \5 k + d(x,y'). 

and if the above estimate holds with x and y interchanged, then the above almost-orthogonality 
estimate (14.4K still holds, but with y replaced by some 7 / G (0,?/). 

Assuming the claim for the moment, we obtain that 


sup \D k/ (f)(z)\ 

^Q k 'r 


< E /*«*.) (/• M 

k £Z V PWal 


^0 


^\k—k'\ri 


X 


fi(kAk') 


+ R„. a «(x‘:+ k ) + V(x*,x*+ K ) V<“»'> + d(x*,x*+ N ) 

As a consequence, we have 

1/2 


J2 su p \ D «(f)( z )\ 2 XqK'+n (x) 
k' al ^y+N zeQ%+ N 

<c{E E E E "W‘)|( f’ 

t /5v AT l-<{Z. r 77 (7)/]<•. 


4>a 


k' a'£%' k '+ N k£Z ae^ k 

1 


\fiAQ 


X 


fi\k-k’\ri 

fi(kAk') 


Xqk'+N (a:) 


1/2 


V 5 (k'Ak)(x k ) + Vtffc/Ak)(a) + v(x k ,x) V 5 ( feAfc ') + d(x k ,x) 

Using the same estimate as in ED, pp.147 148 (sgg cilso Lghiiiici 2.12 111 |HLL2 j j 5 wg obtciiii 

1 / k') 


E #*(<£) 




V S (k'A k)(x k ) + V S (k'Ak)(x) + V(x k ,x ) Vh( fcAfc ') + d(x k ,x 


7 / ib 

f, ^ 




fc]oj(l— 1/r) 


{^(E|(/. 








hr 


1/r 
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where M. is the Hardy-Littlewood maximal function on X and < r < p. 

Thus, by the Fefferman-Stein vector-valued maximal function inequality with p/r > 1 
(see ESI), we obtain 

1/2 

LP(X) 

,21 X /2 
Mr l 


{5Z SU P^ \ D k'(f)(z)\x Q v+N(-)} 

k' o 

< C 


v a'e ar k '+ N z ^Q k a ' +N 


E E lwt/>&*(-)fj 

k ae& k 


LP(X) 


It remains to show the claimed almost-orthogonality estimate (14.41) . We first consider the 
case k > k'. Applying the cancellation property for ijj k (x) yields 

/ V £(0 







f *l>a( x ) r 


JX A JKQa) 


< 


\^*( x )\ 


D k >(x,z) - D k i(x k a ,z) dp(x) 
D k >(x,z) - D k >(x k ,z)\ dfj,(x) 

D k >(x,z )| + \D k >(x k ,z)\ dfi(x) 


JWi a /KQaJ 
+ f \^a( X )\ 

JWi V ll(Q k 

=■■ Ui + U 2 , 

where W\ := (i 6 I : d(x,x^) < {2A Q )~ 1 {8 k ' + d{x^,z))} and W 2 X \ W\. 

Similarly to the estimate of (A) : in the proof of Lemma [3.61 for Ui, using the size condi¬ 
tion (Definition 13.11(4 )) on ^ k (x)/yjp(Q k ) and the smoothness condition (Lemma l3.61( ii) ) on 
D k '(x,y), we obtain that for all z, x^, +N E Q^i +N , 

5 k 


Ui < C 


x 


Wi Vs k ( x a) + V ( X a, X ) ^ + d( X a, X ) 
d(x,x k ) Y 1 


5 k ' + d(x k , x" +N )) V s y (x k ) + V(x, x k ', +N ) \5 k ' + d(x k a , x*+ N ) 


dp(x) 


< C8^ k ~ k '^ 11 


8 k 


'w 1 V S k(x k ) + V(x k ,x)\S k + d(x k ,x) 

1 


r -v 


dp{x) 


x 


V s y (xi) + VUM) + V(xi,x* +N ) W' + d(x*,x* +N )- 
for T > r) and 7 > 0 . 

The estimate for U 2 is similar to the proof for (A ) 2 as in Lemma [3761 Specifically, we have 
U 2 <C 


[ 1 

f 8 k 

Y 

1 1 

( ^ ) 

Jwi V sk (x*) + V{ x %, x ) 

+ d(x k , x) 

) 

W sl/ (z) + V(x,i*, +N ) 

+ d(x,x%f +N )' 


+ 


V (a£) + V ( x k a , x*+ N ) \8* + d{x k , x k „, +N ) 


< C8^ k ~ k '' >v 


V s «(xi) + Vplx*?") + V(xlx*+ N )\g* + d(xl,x« +N 


dp(x) 

5 k ’ 


These estimates of Ui and U 2 establish the claimed almost-orthogonality estimate (14. 4 p 
when k > k 1 . The proof for the case k < k’ is similar. This completes the proof of the 
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almost-orthogonality estimate (14. 4 j) . and hence the proof of the first Plancherel-Polya in¬ 
equality (14.2)1 . 

To show the second Plancherel-Polya inequality (14.311 . we need the following result about 
the operator T x , as mentioned in the outline of the proof of Theorem 14.41 


Lemma 4.6. (Properties of T N ) Suppose that f £ L 2 (X ) and < p < oo, where to is the 
upper dimension of (X,d,/T). Let N be a positive integer. In each cube Q k c f rN , fix a point 
x^ +N ■ Define the operator T N by 

(4.6) T N (f)(x):=Y, E ^Ql +N )Ddx,x k * N )D k U)(^ N ). 

k aeJTHJV 

Then the following assertions hold. 

(i) T n is bounded on L 2 (X). 

(ii) There exists a constant C independent of f and of the choice of x^ +N such that 

\\S(T N (f))\\ LP{x) <C\\S(f)\\ LP{x) , 


where S is the discrete Littlewood-Paley square function as in Definition f.l 


(iii) If N is chosen sufficiently large, then T x is invertible and there is a constant C 
independent of f and of the choice of x^ N such that 


(4.7) 


l|S(7y(/))lli*(x) < C||5(/)||i» (x) . 


We defer the proof of this technical lemma until after the end of the proof of Theorem 14.41 
We now continue the proof of the second Plancherel-Polya inequality (14.311 . Choose N 
sufficiently large that T N is invertible and (14.71) holds. For / £ L 2 (X ), write / = T^T N f. 
Applying Lemma 14.61 we find that 

,2fV2 ^ _ , ., 2 W/ 2 


E E \(€,) 

k a£^ k 


Lp( X) 


(4.8) 


< C 


E E \X a ,T^T N f)x Qi (-)\ 

k 

E E M’ T ‘ v/>x<j*(o 

k oey* 


LP{X) 


1/2 


Lp(X) 


By the definition of T x (f), we have 
/ fpa 

Vh(QaJ k' a '&3X k '+ N 

Therefore, for each fixed rf £ (O,?/), 

12 


Xn!) = E E P(Q e X){^yO t '(-,^ +N ))D t '(f)(x k+N ). 


\(^ T Nf)XQ k (x)\ 

k a G9 k 


11/2 


<p{EE[E E uqW) v-*' 1 ”' 

k ae^ k fc'SZ a'& 3 C k '+ N 
§(k/\k') 


v«)+v(4 + ") + n<,=% +N ) 


.4<* At ') + <i(iJ,i/, +A ') 

s c {EE[E <M "' E /*««' 


D k'(f)(Xa' +N )XQ k (x) 
+JV) 


211/2 


V s t'(x) + V S k'(x*, +N ) + V(x,x*, + ") 


t J +JV\ 


k a.G& k k'GZ 


a’£% k '+ N 
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X 


£(kAk') v 'y 

S^) + d{x,x^ +N )) Dk '( f ^ x ° ,+N ^( x \ 


1 2 >, 1/2 


By the same estimate from EH as in the proof of (j4.2f) above, we have 

1 


E m<?W) 

a'e%~ k '+ N 


X 


Vi-'W + v si ,(x*+ N ) + V(z,x* +K ) 

Sj{kAk') N 7 

^'(/)(4 +JV ) 


^(fcAfe') _|_ x « ( 


fc'+JVi 
a' ) 


^/g^fc'+JV 


1/r 


where A"! is the Hardy-Littlcwood maximal function on A" and < r < p. Note that the 
above inequality still holds when the point x h a xN on the right-hand side is replaced by an 
arbitrary point 0 in Q k a + N , and therefore also holds when the expression | D k i(f){x k a f rN )\ r on 
the right-hand side is replaced by the inhmum of \D k \z)\ r over all z G Q^/ +N ■ Thus, we have 

KV£,/)Xq* 7 


x 


k Q ef fc 


< fi\k-k'\efi[(kAk')-k\u}(l-l/r) 

^ k a£^ k '-k'GZ 


X { M { Yl ^ v \ D k'(f)( Z )\ r X Q k’+N(-))(x) 

V \ “I - -'’ ^ Ot! / 


yg^-fc'+iV 


1/r 


XQ k { X ) 


1/2 


Applying the Fefferman-Stein vector-valued maximal function inequality with p/r > 1, 
from [FSj . we obtain 


\^tf)xQ k a { x )\ 

k ae^ k 


1/2 


LP(X) 


< C 


Y 1 l< ^a’ T Nf)XQk(z)\ 

k a^ k 


1/2 


LP( X) 


< C 


< c 


fc'ez 


- rcn^-i 

'e& k ' +N 


{E E EU Di ' (/)W 

k’ a ' £ %-k>+N z ^ a , 


1/2 


,\ / xl 

2/r-, 1/2 

))o*o 

} 

LP(X)’ 



LP(X) 


which implies that the second Plancherel-Polya inequality (14.311 holds for / G L 2 (X). The 
proof of Theorem 14.41 is complete, except for the proof of Lemma [4.61 □ 


It remains to prove the technical lemma used in the preceding proof. 

Proof of Lemma (i) Fix (VgN. We show that the operator T N is bounded on L 2 (A). 
Write T N (f)(x) = E k (f)(x), where the kernel E k (x,y ) of E k is given by 

E t (x,y):= E p(Q k a N )Ddx,x k * N )D t (x k p N ,y). 
aesr k + N 

This kernel E k (x,y) satisfies the same decay and smoothness estimates (13.91) and (13.101) as 
D k (x,y) does, with bounds independent of x^ +N , as can be shown by a proof similar to that 
for D k (x, y). Moreover, SxM x , y) clp(y) = 0 for each x G X and f x E k (x, y) dpfx) = 0 for 
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each y £ A". Therefore the Cotlar-Stein lemma can be applied to show that Ty is bounded 
on L 2 (A). 

(ii) Suppose that / £ L 2 (X) and < p < oo. Then by the definition of D k and the wavelet 
reproducing formula (12.131) . we have 

f(x) = J2 D t D t (f){x) = J2 E M k * N )D k {x,x k * N )D k (f)(x k + N ) 

k k aefHJV 

k k a&3(' k + N 

= : T N (f)(x ) + R N (f)(x), 

where are arbitrary hxed points in Qq + jv . 

Since T N = I — R N by definition, to show (ii) in Lemma 14.61 it suffices to show that 

(4.9) ||S(fl w (/))lk»m < Ci’"'||5(/)|| 1 , m . 

For then 


l|S(r„(/))||„ ( . Y) < ||S(/)|L, m + ||S(fl K (/))||„ m < (l + CV' v )||5(/)||„ m , 


as required. 

To establish (14. 9p . we write 

Ajv(/)W = E E / {D t (x 1 z)D i (f)(z)-D k (x 1 x k + N )D k (f)(x k + s )}dv(z). 

X- J ' J n k+N 

k a&ST k + N Qot 

Thus the kernel R^^x^y) of Rn is given by 


Rn(x,v) ■■= ! [D k {x,z)D k (z,y) - D k (x,x^ +N )D k (x^ +N ,y)\dfx(z) 

LL — n k + N 

k a&3C' k + N 

= X] [ l D k{x,z) - D k (x,x k a +N )}D k (z,y)dy(z) 

z — n k + N 

k a&3C' k + N 

+ Y1 f D k(x,x k a +N )[D k (z,y) - D k (x k a +N ,y)]dp(z) 

LL — n k + N 

k ae^r k + N 

= : R$(x,y) + R%\x,y). 

Note that by the same proof as for 'Ey, both R$ and R$ are bounded on L 2 (A"), and 
therefore the inner products (ip k , R^\f)) and (ip k , R^\f)) are well deffiied. To estimate 
I|S(< ) (/))IIW), we write 


\\S(Ri\m\L P{ x) 


{EDAsh/i^H l 2 } 1/2 

k a€?y k 


LP(X) 


By the L 2 (A")-boundedness of R^ and the wavelet reproducing formula (j2.13[) for / £ L 2 (A), 
we have 


/ a 



IX 


1>a(x) 

VKQa. 


Z^ z^ 

k 1 a'e^ k '+ N 


Q 


k’ + N 


[D k f(x,z) - D k , 


(x,x k :+ N )} 
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// // 


x/3'M£ E (V£", (y) d M d K z ) d K x ) 

k"ez a " e ^k" 


' x 


Y Y v(q 

k" a "e& k " 

Y Y mq 




./>£ 



^ai x ) rT t \ ^"(y) j , \ j , \ 

Dk'{x, y) dn(y) dn(x) 


»(Q a ») * JXJX 

' rv 


V^iQa) 


M k n ") 


A/ ^(0 F 

J/\ / / .—\ 7 - \ l P l\ /I 


'e& k 




where 


and 


D kl (x,y) := V' / [ZA(a;, z) - A fc /(a:, ^', +7V )] D k >{z,y) dfi(z) 

L ' //ofc' + iV 

a'e«r fe, + JV 


: = 


■=Y D k\x,y) 

k ’ Jx 


^(y) 


y(Q 


d M- 


We now show that yPa(-) /\/ lAQ k a )i A" (’) A") / satishes an almost-orthogonality estimate 

similar to (14.411 . by following the philosophy of Remark [4.51 Recall that ^a(') / V d(Qa) is a 

test function. It remains to show that the function F k » (x, x k „) satishes a size condition, a 

Holder regularity condition, and cancellation. 

// 

Next, it seems unlikely that F k "(x,x k „) satishes the Holder regularity condition (I3.10p . 

However, as noted in Remark 14.51 it suffices to establish the weaker Holder regularity condi- 

_ ^ ,// 

tion (14.51) . which we now do. To begin, we show that F k "(x,x k „) satishes 


(a)' \F h .,(x,x\,)\ <CS” N 


V, h "(x) + V(x,xK',) + d(x,xK) 


S k 


for all 7 G (0, r/), and 

(by 


<CS r, N (d^)Y 


+ 


x 


\y 5 u"{ x ) +v{x, xK) +d(x,xK) 


5 k 


" \ V // 


V 6k " (x') + V (x', x™„) + d(x', x%) 


" \ \ rl" 


for all ?/ G (0,?/). 

To prove (a) and (b) , we hrst show that D k i(x,y) satishes the same estimates (13.9j) 
and (13.101) as D k (x,y), but with the constant C replaced by CS 7 ' N , that is, 


(a) \D k i(x,y)\ < C5 vN - 


( 6k ' \ 

7 

l ; 

i / A' x 

x, y ) \5 k ' + d(a 
/ d(x,x') \ 

r ',y)' 

V 

\<W + d(x,y )/ 

V s y(x) + V(x,y) \S k ' + d(x,y)J 


7 
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for d(x, x') < ( 2A 0 ) 1 maxl^’ + d(x, y ), 5 k ' + d(x ', y )}, and 
(c) \D kl (x,y)-D kl (x, y ')\<CS^ N \ 


d(y, y') w 


-S k ' + d(x, y)) V S k> (x) + V(x,y)\8 k ' + d(x,y) 
for d(y, y') < (2A 0 ) _ 1 maxl^ + d(x, y), 8 k ' + d{x 1 y')}. Moreover, 


D k '(x,y) d/i(x) = 0 and / D k >(x, y) dy(y) = 0, 


lx 


ix 


for all y E X and all x G X, respectively. Indeed, note that [D k /(x,z) — D k >(x, ^ ,+iV )] 
satisfies the same estimates (13. 9 p and fl3.10p as D k /(x,z) does, but with the constant C 
replaced by CS vN . Therefore, the proofs for (a), (b) and (c) follow from a similar proof to 
that for Lemma 13.61 As a consequence, the almost-orthogonality estimate (1 1.1 li holds for 

D k f(x, ■),ip k "(-)/\/li(Q k ii) V We omit the details. 


Now, to verify the estimate in (a), applying this almost-orthogonality estimate yields that 

rA-) 




k’ 


d(Qt") 


< 


< 


k \r] 


§ (k'Ak ) 


cs^Ysw- 

V v s o b'Afc") 0) + y (*, xk a ") V + d{x, xK ) 


CS vN (5l fc - fc, l(7-7)___ 

v V 8k" ( x ) + V ( x i xk J >) ^ S k " +d(x, xK ,) 


< C5 vN 


V s t" 0 ) + V ( x , x a") ' d k " + d(x, x k J ,) - 


for all 7 G (0, rj). 

Next we show the estimate in (b) 7 . Note that 

\D k ix,y)-D k ,(x',y)\ < 'X-)” 


X 


5 k 


+ 


V 5 v{x) + V(x,y) \5 k ' + d(x, y)) V 5 v(x') + V(x,y) \5 k ' + d{x’,y) 


71 


and 


ix 


D k ’(x, y ) - D k i(x', y) dy(y) = 0 . 

Therefore, as pointed out in Remark 14.51 we obtain for k! > k" that 
E k”( x , xk „") ~ E k"(, x 'i x l") 


< Y \([ D k'( X ’ 0 ~ D k'( x ', •)]> 


ib k 

' rt 


k’ 


< C5 vN E fi\k'-k" 1 1) 

1 


v M k n ») 

-k"\r,f d ( X ’ X ')V' * 1 


k' 


5 k ' 


L y g k" ( x ) + V(x, x k „ )^S k + d(x, x k „) 


" \ \ c-l." 


+ 


6 k 


V k " (x') + V (. x ', x k "„) A S k " + d(x', x k "„) 


71 


5 k 
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< C5 vN <*'*'" 

k’ 

+ 


-k"\(r t -v')( d (x,x') y' 


Sk " ' L Kfc")0) + V ( x ’ xk ") V Sk " + d( K X i X t") 


a 


V S k")(x') + V(x',x k „) ' 5 k + d(x',xK) 


>> 
V 


" \ V " 


/ Z/' 

a 


< 


+ 


/d(x, x') 

r' 

r 1 

( ^ ) 

V 5 k " 

) 

( X ) + V ( X ’ X a") 

^8 k ” + d(x,xK) d 


1 

n 

( 8 k 

VI 



K fc « (a/) + V(x, x k "„) ^ 8 k " + d(x', x k "„) 


5 k 


for all rj' G (0, 77 ). 

For k 1 < k" , we have 


< y' / i[-Dfc'(® > j/)--Dfc'( a;, >j/)]i 
fc '< fc " 

J d(;r, a;') W r 1 

— ' 


^ k >> ( y ) 


< 






dKy) 

V 


+ 


V s y (x) + V{x, y) \5 k ' + d(x, y) 
S k ' \ 71 1 


W) + V(x',y) \8 k ' + d(x',y)J J (y) + V(y,xK) + d(y,xK,) 


dy(y) 


< 


C§ vN( d( X i X ') y jr I(k"-k') ' 

8 ....// 




a-' 


+ 


V*' (x) + V(a, a£„) v <5 fe ' + <Z(x, ar „) 

V \ 7’ 


< 


V VO + W V„) ^ 8 k ' + d(V, V„) 

C5 vn 58 i-i)(k-y) ■ 




V*" V) + v(x,x k „) + d(x,x k „) 


S 


k" 


+ 


5 


fc" 


< 


V^fc"(x') + V{x',x\) ^5 k " + d{x',x k a „) 
d(x,x')\v[ 1 


8 k ' '-V 5 k"(x) + V(x,x k „) + d(x,x k „) 


zk" 


s :k" 


V 5 k"(x') + V(x',xK) + d{x',xK) 
for all 7 G (0, 7 ). 

With the almost-orthogonality estimate in hand, the same argument as for (14.21) . via the 
estimate from (ED and the Fefferman-Stein vector-valued maximal function, yields 


= ||{E E KC^t/MoitOf} 


1/2 


k ote<3f k 


LP(X) 
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< CS vN 


E E K<.(/)>x«*( 

k 


1/2 


LP (X) 


= CT7|S(/)||„„>. 


A similar proof shows that ||*S'(i?^' ) (/))||i / p(x) < C5 7lN \\S(f)\\LP(x)- Therefore (14.9p holds: 

\\S(R N (f))\\ LP{x) <C5* N \\S(f)\\ LP{xh 


as required. 

(iii) Consider the Neumann series (T/v) _1 = (/ — R N ) _1 = By (14.911 we have 

OO 

l|S((r w )- 1 (/))||„, A . ) < E IIS((flU i (/))ll«.(x ) < (i - ct’' n )- 1 ||S(/)||„ (a) , 

i =0 

as required, if N is chosen sufficiently large that C5 vN < 1. 

This completes the proof of Lemma 14.61 □ 

We turn to the product setting. 


4.2. Product square functions via wavelets, and Plancherel—Polya inequalities. 

We now assume that X = X\ x X 2 where each X t is a space of homogeneous type as above. 
In this subsection, (rci, rc 2 ) denotes an element of X 1 x X 2 . 


Definition 4.7. (Product square functions) Take fa G (0,r/j) and 7 * > 0, for i = 1, 2, and 
consider / G (G(fa, fa] 7u 12 ))'- 

(a) The discrete product Littlewood-Paley square function S{f ) in terms of wavelet coefficients 
is defined by 


(4.10) S(f)(x i,x 2 ) := 


{E E E E f)x Q ^Mx Q ^M 

fel fc 2 a 2 G^ k 2 



where V ; a)' 0 a 2 = ® V’aa with acting on the W variable for z = 1 , 2 , and (a;*) := 

X Q k d x i)fM(QaJ~ 1/2 - 

(b) The continuous product Littlewood-Paley square function S c (f ) in terms of wavelet oper¬ 
ators is defined by 


(4.11) 


S c (f)(x i,i 2 ) := 


EE 

fci fc 2 


D kl D k2 (f)(x 1 ,x 2 ) 


A 1/2 

/ ’ 


where D ki := J2 ai ^ k i V’a- for * = 1, 2 , and D kl D k2 := D fc] ® D fc2 . 


The main results of this subsection are the following product versions of the Littlewood- 
Paley theory and the Plancherel-Polya inequalities. 


Theorem 4.8. (Product Littlewood-Paley theory) Suppose fa G (0, 77 *), 7 * > 0, and 
max{} < p < 00 , where u>i is the upper dimension of X l7 for i = 1,2. For 

all f G (G(fa, fa] 71 , 72 )) , we have 

\\SU)\\lp(X 1 xX 2 ) ~ \\S c (f)\\LP(X 1 xX 2 )- 

Moreover, if 1 < p < 00 , then 

\\S(f)\\ L P(X 1 xX 2 ) ~ \\S c (f)\\LP(X 1 xX 2 ) ~ ||/||l, p (V 1 xX 2 )- 
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Theorem 4.9. (Product Plancherel-Polya inequalities) Suppose f G (0 ,r}f), 7 \ > 0, and 
max{ M2+m } < P < 00 > where u>i is the upper dimension of X i; for i = 1,2. Take N { , 

N 2 G N. Then there is a positive constant C such that for all f G (G(Pi, fd 2 ] 7 i, 72)) , we have 


E E E E 

k l a' 1 eJ' fc l +iv l k 2 a' 2 &X k 2 +N ‘2 

sup 

, , „kt +N-I „k'+Nn 

( 21 , 22 ) 6*9 ) 1 xQ 2 , 


\Dk[Dk' 2 (f)(zi, z 2 )\x fc'+JVi (Ox k ' 2 + N 2 (•) 

Q 1 a 2 


1/2 


LP(X ixJfj) 


( 4 -i2) <c-|{e; e e e i<«./)x oSi (-)x oa (or} 


1/2 


fel aiS^ fc l fc 2 o 2 6^ fc 2 


Lp(A'ixX 2 ) 


Further, suppose JVj and N 2 are sufficiently large positive integers, to be determined during 
the proof below. Then there is a positive constant C such that for all f G (G(/3i, /3 2 ; 71 , 72 )) , ? 
we have 


E E E E |(V£fe/>x Q *i(-)x 

fci a 1 G^ fc l k 2 ot 2 g @^ fc 2 


1/2 


/ 


Lp(X 1 xX 2 ) 


< C 


E E E E 

fc, i a' 1 e^r fc 'i +JVl fc 2 a^e^' fc 2 +JV 2 


(4.13) 


inf 


, , „k( + iVi „fci + JV 2 

( 2 l, 2 2 )GQj xQj 


|-Dfc'-Dfe'(/)(^ij 22 )) X Q fc'i+JVl (■)X r) fci+^2 


QJ 


1/2 


Lp(A7xA 2 ) 


Proofs of Theorems \fP\ and \f£\ The proofs of Theorems 14.81 and 14.91 are analogous to those 
for the case of one factor. As mentioned in that case, the proofs of Theorems 14.31 and 14.41 
follow from the almost-orthogonality estimates, namely the claim (14.41) . To see that these 


proofs can be carried over to the product case, we make an observation analogous to the 
claim (14.41) . as follows: 


D 


k 1 D k2 (x 1 ,x 2 ,-,-),'if* 1 l f'£(-,-,y 1 ,y 2 )^ = {D k fx 1 ,-),fl\(-,y l )){D k2 (x 2 ,-),f k f 2 (-,y 2 )), 


which together with the almost-orthogonality estimates for the one-factor case yields the 
desired almost-orthogonality estimates for the product case. For the product case, all esti¬ 
mates analogous to those in (a)-(c), (a) -(b) 7 , (j4.2[) and (j4.3[) follow similarly. We omit the 
details. □ 


5. Product H p , CMO p , BMO and VMO, and duality 

In this section we define the Hardy spaces II P . the Carleson measure spaces CMO p (includ¬ 
ing the bounded mean oscillation space BMO = CMO 1 ), and the vanishing mean oscillation 
space VMO, in the setting of product spaces of homogeneous type. Both H p and CMO p are 
defined here for p in the range max{^-, 777 ) 7 } < p < 1 , where u 7 is the upper dimension 
of Xi, for i — 1, 2. We prove that CMO p is the dual of H p , and in particular that BMO is 
the dual of H 1 , and also that H 1 is the dual of VMO. 

We develop this theory in the product case with two parameters. The generalization to 
k parameters, fc G N, is similar to the two-parameter case, while the specialization to one 
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parameter is immediate. Note the difference from the Littlewood-Paley theory developed in 
Section [4] above; there it was necessary to develop the one-parameter theory first, then to 
pass to the product case by iteration. 

Fix (3i G (0, rji ) and 7 * > 0, for i = 1,2. 

OO o 

For brevity, we denote by G and ( G)' the test function space G{fi 1 , ^ 2 ; 7i, 72 ) and the space 
of distributions (G(/3i, /3 2 ] 7i, 72 )) , respectively. 

We are now ready to introduce the Hardy spaces H p ( Ad x X 2 ) and the Carleson measure 
spaces CMO p (X 1 x Ad). In this section, (x\,x 2 ) denotes an element of X\ x X 2 . 


Definition 5.1. (Hardy spaces) Suppose max{, 7 ^ 7 ^} < P < 1, where 07 is the upper 
dimension of A"* for i — 1,2. The Hardy spaces H p (Xi x X 2 ) are dehned by 

H p (X 1 X X 2 ) := {/ G (GY : S(f) G X X 2 )}, 

where S(f) is the discrete product Littlewood-Paley square function as in Definition 14.71 
For / G H p (X 1 x X 2 ), we define \\f\\ H p(x 1 xx 2 ) ■= \\S(f)\\ L p(x 1 xx 2 )- 


For completeness, we note that as in the classical case, for 1 < p < 00 the Hardy 
space H p (Xi x X 2 ) of Definition 15.11 coincides with L p (Ad x X 2 ). 

We point out that G and hence H P (X\ x X 2 ) D L 2 (X 1 x X 2 ) are dense in H p (X\ x X 2 ). 
Indeed, if / G H p (Xi x Ad), then by Theorem 13.41 the functions 


fn{xi,X 2 ):= 

\k 1 \,\k 2 \<n aie & k i 02G^ fe 2 


belong to G■ Moreover, 
S(f ~ f n )(xi,x 2 ) < 


|(V’aiV’a 2 2 ,/)X Q M(^l)XQ^(x 2 ) 

|fcl |>n or \k 2 \>n a 1 £^ k l a 2 e ^ k 2 


2 11/2 


Therefore ||5(/ 
H p ( Ad x X 2 ). 


o 

fn ) 11 up (Yi x x 2 ) tends to zero as n tends to infinity. Hence G is dense in 


Definition 5.2. (Carleson measure spaces, and bounded mean oscillation) Suppose that 
max{ i ^+^ 2 } < P < 1) where c <7 is the upper dimension of Ad for i — 1, 2. We dehne 
the Carleson measure spaces CMO p in terms of wavelet coefficients by 

CMO p (A 1 x X 2 ) := {/ G (£)' : C p (f) < L°°}, 


with the quantity C p (f ) dehned as follows: 


(5.1) 


iy7.i := sup 
n 



E 

fcl ,fc2 ez,ai 60^1 ,a 2 e& k 2 



where D runs over all open sets in Ad x X 2 with hnite measure. 

The space BMO of functions of bounded mean oscillation is dehned by 


BMO(Ad x X 2 ) := CM 0 1 (A 1 x X 2 ). 


The main result in this section is the following. 
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Theorem 5.3. Suppose max{} < p < 1, where uii is the upper dimension of Xi 
for i — 1, 2. Then the Carleson measure space CA10 p (Ai x X 2 ) is the dual of the Hardy 
space H p (X 1 x X 2 ): 

(H p (X i xl 2 ))' = CMO p (Ad xX 2 ). 

In particular, 

(H\X i x X 2 ))' = BMO(Xi xX 2 ). 

To prove Theorem 15.31 we follow the approach developed in (HLL2] ; see also [HLLlj . We 
first recall the definitions of the product sequence spaces s p and c p for 0 < p < 1. These 
sequence spaces are discrete analogues of H p (Xi x X 2 ) and CMO p (Ad x X 2 ) respectively. 
The space s p is defined to be the set of sequences s = {sp}p of real numbers such that 

(5-2) NI* := \\{Y1 1 p(R)~ 1/2 srXr(-, •)| 2 } 1/2 || L p (Xi xx 2 ) < 00 > 

R 

where R runs over all dyadic rectangles in X l x X 2 . The space c p is defined to be the set of 
sequences t = {tp}p of real numbers such that 

(5.3) \\t\lcp := sup ( N 2 ) ' < 00 , 

n V(n)" 1 Kta J 

where runs over all open sets in Ad x X 2 with finite measure, and R runs over all dyadic 
rectangles contained in 

We emphasize that in the above definitions of s p and c p , the expression “all dyadic rect¬ 
angles if” indicates the rectangles of the form R = Q k c f x Q k ff 2 for all ki € Z and cq G 3H ki 
for i — 1, 2. 

The main result about the sequence spaces s p and c p is the following duality result. 
Proposition 5.4 ([ HLL2 ]). For 0 < p < 1, (s p )' = c p . 

We now introduce the lifting and projection operators Tl and T p , as follows. 

O 

Definition 5.5. For / e ( G)', the lifting operator Tl is defined by 

(5-4) {(TlJ)r}r~ {W'SV>S./)} s . 

where R — X Q k ^ , k 2 £ Z, Or G ty kl , a 2 € l ty k ' 1 are dyadic rectangles in Ad x X 2 . 

Definition 5.6. Given a sequence A = {A/;} of real numbers, we define the associated 
projection operator Tp by 

(5.5) T P (X)(xi, x 2 ) := ^R-^ k a \{ x M k a \{x 2 ). 

R=Q k a \xQ k J 2 , 

fci to ez,ai ,a 2 er fc 2 

From the definitions of the lifting and projection operators Tl and Tp, it follows that 

o o 

f = Tp o T L (f) in the sense of the test function space G and of the distributions (G)'■ That 

O 

is, Tp o Tl is an identity operator on the distributions (G)'- 

Next we give two auxiliary results which will be used in establishing the duality in Theo¬ 
rem 15.31 
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Proposition 5.7. Suppose max {■ ^ 2 +m } < P — where uii is the upper dimension 
of Xj for i = 1,2. Then for all f G H P (X 1 x X 2 ), we have 

( 5 -6) \\T L (f)\\ s p < \\f\\Hp(x 1 xx 2 )- 

In the other direction, for each s G s p we have 

( 5 . 7 ) \\Tp(s)\\hp(x 1 xx 2 ) < ll s ll^- 


Proof. Inequality (15.61) follows directly from the definitions of H p (Xi x X 2 ) (Definition 15. lj) 
and the sequence space s p (formula (15.21) 1. 

We now prove (15.7j) . For each s G s p , by the definitions of H p (Xi x X 2 ) and Tp(s), we 
have 


\\Tp(s)\\ H p( Xl xX 2 ) - II*S'(^ 1 p('S))||lp(ATxA 2 ) 


E 


R=Q7 1 xQ7 2 , 


E 

k 1 

fciV^GZ.a'G^'l ,a' 2 e^ k 2 


k f k r 

R'=Q ) xQ 2 






Xr> fc i (^l)Xr) fc 2 (^ 2 ) 

W Q ~ 


2 -> 1/2 


E 


Lp(XixX 2 ) 

■Sr • Xn k 1 (zi)£y=2 (*2)l^ 

Wen x Wq:2 J 


'n 1/2 


R=Q,fxQ7 2 , 

fcl,A:2GZ,aiG^ fc l ,a 2 G^ fc 2 


LP(ATxA' 2 ) 


— ll S l|s p ) 

where the third equality follows from the orthogonality of the bases {f^\ } and }. □ 

} < p < 1, where cuj is t/ie upper dimension 


Ul 


UJ 2 


Proposition 5.8. Suppose max { 

r ^ 1 £Ji+r;i ’ w 2 +?7 2 

of Xj for i = 1,2. For all f G CMO p (Ab x X 2 ), we have 
(5-8) \\T L (f)\\ c P < C p (f). 

In the other direction, for each t G c p , 

( 5 . 9 ) C p (Tp(t)) < ||<||„. 

Proof. Inequality (15.81) follows directly from the definitions of CMO p (A r ! xX 2 ) (Definition [5T2]) 
and c p (formula (15.3f) ). 

We now prove (15.9ft . For each t G c p we have 


C p (Tp(t )) = sup 


l p(D) 


E 


TpM>| 2 } 


1/2 


/ 


R=Q 7 1 xQ a 2 2 cn, 

fcl,fc 2 6Z,aiGr fe l ,a 2 G@^ fc 2 


= sup 
n 


l fi(Q) 


E 


HR 


1/2 


R=Qa 1 1 X < 3-2 Cn > 

fcl,fc 2 SZ,oiG#' fe l ,a 2 GS^ fc 2 


< 


| CP) 


where the second equality follows from the orthogonality of the bases {V'm} and {f’al}- ^ 
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We would like to point out that thanks to the orthogonality of the wavelet basis from [AHi, 
the proofs given here of (15.71) and (I5.9j) are much simpler than those given in |HLL2 j. 

We are ready to prove the duality [H p {X\ x X 2 ))' = CMO p (Ah x X 2 ). 

Proof of Theorem \5.,‘A Suppose max{ , ^ 2 +m } < P — 1- We s h° w that there exists 
a positive constant C such that for each g E CMO p (Ah x X 2 ), 

(5-10) \(f, g) \ < C\\f\\ H p(Xixx 2 )C P (g) 

for all / E G- It follows that CMO p (Ah x X 2 ) C (H P (X i x X 2 ))\ since G is dense in 
H p (X { x X 2 ). 

o 

To prove inequality (15.10b for each f E G and g E CMO p (Ah x X 2 ), by the reproducing 
formula (13.261) we have 

(f,g) = (V’ai/KV’alife d) 

fci,fc 2 GZ,aiG^ fcl ,o;2G^ fc2 

= £( t M))r-(Tl(</))r- 

R 

where Tf,(/) and Tjfg) are the lifting operators as in Definition 15.51 
Then, by Propositions 15.71 and 15.81 we obtain 

K/.J)I < \(Tdf),T L m < C\\f\\ H p(x 1 xX2)Cp(g)- 

Conversely, suppose l E (H p (X 2 x X 2 )) . Let l\ := l o T P . By Proposition 15.71 we see that 
h E (s p )', since for each s E s p , \h(s)\ = |/(T P (s))| < C\\l\\ ||Tp(s)|| ffP(XlxX2) < C\\l\\ ||s|| sP . 
Now we have 

l{g)=loT P oT L (g) = h{T L (g)) 

o 

for each g E G- So by Proposition 15.41 there exists t E c p such that h(s) = (t, s) for all s E s p 
and ||t||cP ~ ||/i|| < ||/||. Hence 

l(g) = {t,T L (g)) = (T P (t),g). 

By Definition 15.21 and Proposition 15.81 we obtain that ||Tp(t)||cMo p (Xixx 2 ) ||i|| c p ||(||. 

Hence (H P (X 1 x X 2 ))' C CMO p (Xl x X 2 ). □ 

Now we introduce the space of functions of vanishing mean oscillation. 


Definition 5.9. (Vanishing mean oscillation) We define the space VMO(Ah x X 2 ) of functions 
of vanishing mean oscillation to be the subspace of BMO(Ah x X 2 ) consisting of those 
/ G BMO(Ah x X 2 ) satisfying the three properties 


(a) lim sup < 

O: ^(Cl)<5 *- 


(b) 

(c) 


r 1 


M n ) fc2 

R=Q,f xQ„|cn 

fcl,fc 2 €Z,a:i63'' i: i ,a 2 €£^ fc 2 




,2-1 1/2 _ 


a\ t ol2 


lim 


sup 


' rCO Cl: diam(S7)>JV ^ h(^) 


E 


= 0 ; 


O'! 1/2 

i<<e/>r} = °; 


R=Q a \xQ a |cQ, 
fci,fc 2 ez,aie®" fcl ,a 2 e^ fc 2 


lim sup 

JV ^ 00 Cl: Qc(B(a:i,A r )xB(x 2 ,V)) c 
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{-4jy E l«e/>| 2 } 1/2 = 0, where 

R=Qa\ xQa 2 2 Cfl 
ki to ez,Qi ,a 2 e^ fc 2 

xi and X 2 are arbitrary fixed points in Ad and Ad, respectively. 

We now show the duality of VMO(Ai x Ad) with H x (X\ x Ad). 

Theorem 5.10. The Hardy space H l (X\ x X 2 ) is the dual of VMO(Ai x X 2 ): 

(VMO(Ai x Ad))' = hT 1 (Ad x A 2 ). 

Proof. The proof of this theorem is similar to the proof of the duality between VMO and H 1 
on Euclidean space given in Section 5 of jLTWj . Following [ LTW] . we only sketch the main 
steps of the proof. First, we use FW to denote the set of finite linear combinations of terms 
of the form {/(/X ■ -0^}, where {yj^f } are wavelets on X t , i = 1,2, as in Theorem 12.21 Second, 
from Definition 15.91 we obtain that VMO (Ad x A 2 ) is the closure of FW in the BMO(X! x Ad) 
norm. 

The inclusion H 1 ( Ad x A 2 ) C (VMO(Ad x X 2 ))' follows from the duality of H l (X\ x Ad) 
with BMO(A! x A 2 ), which was shown in Theorem 15.31 The reverse containment follows 
from the fact that FW is dense in H 1 ( Ad x A 2 ) in terms of the H x (X\ x Ad) norm and from 
the following inequality: for / G FW, 

\\f\\m(x 1 xx 2 ) < C sup |(6,/)|. □ 

6eVMO(AixA 2 ), 

lkllBMO(X 1 xA' 2 )= 1 

6 . Calderon-Zygmund decomposition and interpolation on Hardy spaces 

In this section we provide the Calderon-Zygmund decomposition and prove an interpola¬ 
tion theorem on H p { Ad x A 2 ). Note that H p { Ad x Ad) = L p ( Ad x A 2 ) for 1 < p < 00 . I 11 
this section, (xi,x 2 ) denotes an element of Ad x X 2 . 

Theorem 6.1. Let max{ } < P 2 < 1, where u>i is the upper dimension of Xj for 

i — 1, 2. Suppose p 2 < p < Pi < 00 , a > 0, and f G H p (Xi x Ad). Then we may write 

f(x 1 , x 2 ) = g(x 1 , x 2 ) + b(x i,x 2 ), 

where 

g G H n {X d x X 2 ) and b G H P2 (X! x X 2 ) 

are such that \\g\\ p HP1 ^ XiX x 2 ) ^ Ca pi P ||/||^p(x 1 xx 2 ) an( ^ , \\^\\ P HP 2 (x 1 xx 2 ) — Ca P2 p ||/||^-p(x 1 xa 2 )- 

Here C is an absolute constant. 

Theorem 6.2. Suppose max{} < Vi < Pi < where uj t is the upper dimension 
of Xj for i — 1, 2. Then the following two assertions hold. 

(a) Let T be a linear operator that is bounded from H P2 (X ixX 2 ) to L P2 ( Ad x Ad) and from 
H Pl (Xi x X 2 ) to L Pl (Xi x A" 2 ). Then T is bounded from H P (X 1 x Ad) to L p { Ad x Ad) 
for all p with p 2 < p < pi. 

(b) Suppose T is bounded on H P2 ( Ad x Ad) and on H P1 (X d x Ad). Then T is bounded 
on H p ( Ad x X 2 ) for all p with p 2 < p < pi- 


We first prove Theorem 16.11 
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Proof of Theorem 16.1L Suppose that / G H p (X\ x X 2 ) and a > 0. Let kit := {(xi,x 2 ) G 
Xi x A 2 : S(f)(xi,x 2 ) > a2 e }, where S(f) is the discrete product square function defined 
in (14.101) . 

Let 

n 0 := |f? = Q k a \ x Q k a \ : ki, k 2 G Z, a, G W k \ a 2 G fi(R n fi 0 ) < 
and for £ > 1 


7^ := |i2 = Q k \ X Q* 2 , A*, G Z, a, G ^ fcl , o 2 G 

such that ji{R D f^_i) > ——p,(R) and /i (R D kit) < 7rr~k i {^) 

2 A 0 2A 0 

Applying the wavelet reproducing formula from Theorem 3.11, we have 

f(x 1 ,X 2 ) = EE EE 

fei 1 k 2 a 2 g^ fc 2 

= E E 

l - x R=Ql\xQ k c? 2 e-R e 

+ EE (f^a\^l)^a\( X Ml( X 2 ) 

R=Q k a \xQ k a 2 2 en 0 

=■ K x ,y) + gi x ,y)- 


When pi > 1, the LP{X\ x X 2 ), 1 < p < oo, estimate for the Littlewood-Paley square 
function implies that 


Lpi(X 1 xX 2 ) ^ 


< C 


EE (^a 2 , f)X Q *i ( X l)X Q % (X 2 ) } 


2 -j 1/2 


R=Q k a \xQ k a \&K 0 


LP i (XixX 2 ) 


Next, we estimate ||^||^i(x 1 xx 2 ) when max{^,^} < Pi < 1. We estimate the 
H pi (Xi x X 2 ) norm directly. To this end, using the wavelet coefficients of g, we observe that 

,24 V 2 


H p i(X 1 xX 2 ) < 


E E E E ' ] ( x l)X Q ^{x 2 )\ 2 } 


k 'i a' 1 E& ,k i k i a 2 e ^ k 2 


L p i(XixX 2 ) 


< c 


EE (^K 2 2 > /)Xq*i C&OXq^ O 2 ) } 


2 1/2 


R=Q k a \ xQtl&R-o 


LP i(XixX 2 ) 


Thus for all p! with max{ 


OJl CJ2 

LJl+771 ’ UJ2+ri2 


} < pi < 00 , we have 


H p i (Xi xX 2 ) ^ 


< C 


EE | /)Xq*1 (®l)Xg*2 (^2) 


R=Q^i xQa 2 &TZ 0 


2 11/2 


L p i (Xi xX 2 ) 


Claim 1: We claim that 

S(f)(xi,x 2 ) Pl dp,i{xi) dp 2 (x 2 ) 


1 S(f)(xi,x 2 )<a 

> c 


EE , f)X Q H MXqH {X 2 ) 


R=Q k \ xQ k a 2 2 e^o 


2 ', 1/2 


LP1 (Xi xX 2 ) 
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This implies that 

\hpi(X 1 xX 2 ) E C 


' S(f)(x\,x 2 )<a 


S(f)(xi,x 2 ) pi dni(xi) d/i 2 (x 2 ) 


< Ca pi ~ p 


S(f)(x 1 ,x 2 ) p dfj,i(xi) dp 2 (x 2 ) 


J S(f)(x i,x 2 )<a 

< Ca p '- p \\f\\ p HP{XixX2) . 

To show Claim 1, we choose 0 < q < p\ and q < 2, and observe that 
[ S(f)(x 1 ,x 2 ) pl dfjb^xi) dp 2 (x 2 ) 


1 S(f)(x I,x 2 )<a 


in: 


E E E E ^m^alJ)x Q H( x i)x Q ^ 2 {x 2 ) 


0 k 1 ai ^ k i fci a 2 €^ fc 2 


2-> Pi/2 


/ 


dni(x i) d,p 2 (x 2 ) 


> C 


= C 


> c 


E 


ai^a 2 2 )/)Xgfcl(^l)A 0 ^(T 2 )| 2 ) ' d/Xl (>l) dp 2 (x 2 ) 


0 R=Qa L 1 xQ^|e77 0 


L Qo 


'XixX 2 


'AT xX 2 L 


— 

E /)^i( Xl )E Q fc 2 (^ 2 )xng(^ i ,^ 2 )| 2 f 2 dpifa) dp 2 (x 2 ) 

R=Qa\ xQ«! G 77 o 

E 

R=Q« 1 1 xQ«|e 77 0 


M s{{^a\^alJ) q ^ R ) q Xm%)(x 1 ,X 2 ) 


2 . £ 
<? | 2 


Pi 

9 


d/M(xi) dp 2 (x 2 ) 


> C 


'X 1 xX 2 


2, f 


E <« ,/>XQ*i(® 1 )XQ*2(iC 2 )| j 2 d//i(xi)d// 2 (a:2), 


R=Qa 1 1 xQa%eno 


where in the last inequality we have used the fact that D R) > \n(R) for R e 1Z 0 , and 
thus Xr( x 1 i x 2 ) E 29M s (x Rn ng)«(&!,x 2 ), and in the second to last inequality we have used 
the vector-valued Fefferman-Stein inequality for strong maximal functions: 


E M *(^ 


Lp(X 1 xX 2 ) 


< C 


El/* 


Lp(X\ xX 2 ) 


k= 1 fc=l 

with the exponents r — 2/q > 1 and p = pi/q > 1. Thus the claim follows. 

Let be the enlargement of the set f \ given by := {(xi,x 2 ) G X 1 x X 2 : M s (xn e ) > 

(2^ 0 )- 1 }. 

Claim 2: For p 2 < 1, 

IP2 


E (/.<<)</£(* i)<fe) 

R=Q« 1 1 XQ«| G77« 

Claim 2 implies that 


RP2 (AixA 2 ) 


< C(2 t a)P>(ff £ _ 1 ). 


IP2 


lRP2(ATxA 2 ) — 




— 1; 


< c 


’ S(f)(xi,x 2 )>a 


S(f)(xi,x 2 ) P2 dpi(xi) dpi 2 (x 2 ) 
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< Ca P2 ~ p 


lt x 2 )>a 

< Ca P2 - p \\f\\ p HPiXixX2 y 


S(f)(x 1 ,x 2 ) p dii 1 (x 1 )dn2(x2) 


To show Claim 2, note that 


R=Qa\ xQ a \ene 


< c 


P2 

HP2 (XiXX 2 ) 


f)*Q k a \ ( X ^Q k2 (^2) 


H=Qal xQ%&h 


S') 1/2 


LP2 (XixA 2 ) 


Then 

OO 

E< 2 '“) k ms 


f-ij 


<=i 

> 




S(f) P2 (x 1 ,x 2 ) dni(xi) d/i 2 (x 2 ) 


P2 
T 


,_E E E E (V’a 1 1 V’a 2 2 ,/)X^i(^l)X Q fc|(^2)| ] 2 d^X!) dn 2 (x 2 ) 

'(it-lXVli K kl ai&& -k L fci a 2 &W k 2 


> 


<X 1 xX 2 


E 

R=QaixQa 2 2 e^o 

' kl ^a 2 J)X Q * 1 (®l)X g £ (^2)XQ £ _ A o,(^l,^ 2 )n 2 dfj,i(xi) dfi, 2 {x 2 ) 


N P2 

12 1 T 


> c 


'Xi xX 2 


2-, ^ 


(^mfe/)XQ*i(^i)x o *2 0£ 2 ) } 2 d//i(xi)d/x 2 (a;2), 


R=Oai xQa 2 2 eTZ 0 


where the last inequality follows from the fact that if i? G Ab then i? C f^_i, and therefore 
/i(i?D (0^_!\0£)) > This establishes Claim 2. Hence, the proof of Theorem 16.11 is 

complete. □ 

We end our paper by proving the interpolation theorem on Hardy spaces H P (X 1 x X 2 ). 

Proof of Theorem 1 6. H (a) Suppose that T is bounded from hP 2 (Ab x X 2 ) to L P2 (X! x X 2 ) 
and from H P1 (X 1 x X 2 ) to L pi (Ab x X 2 ). For each given A > 0 and / G H p (Xi x X 2 ), by 
the Calderon-Zygmund decomposition we may write 

f{x 1, x 2 ) = g(x 1, x 2 ) + 6 (xi, x 2 ) 

with 

R p i (A'iX.Y 2 ) — P H/IIrp(AixA' 2 ) an< ^ H^lljrP2(AixA 2 ) — ^11 / II #p(Xi xA 2 ) ■ 


Moreover, we have proved the estimates 


and 


pi < n 

H p i (X\ xX 2 ) — 0 


I P2 

I H p 2 (Ai x A 2 ) — ^ 


S(f)(xi,x 2 )<a 


S(f) pi (xi,x 2 ) d/j,i(xi) d/i 2 (x 2 ) 


S(f)(x !,x 2 )>a 


S(f) P2 (x 1 , x 2 ) dll i(xi) d/i 2 (x 2), 
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which imply that 

II^/IIIp(.YixA' 2 ) = 
< 


< 


< 



a p 1 p({(xi,x 2 ) : \Tf(xi,x 2 )\ > a}) da 
a p ~ x p{{{x i,x 2 ) : \Tg(xi, x 2 )\ > a/2}) da 

POO 

a p - 1 /j,({(x 1 ,x 2 ) : \Tb(x 1 ,x 2 )\ > a/2}) da 


p a- 


p—pi—i 


' S(f)(xi,x 2 )<a 


S(f) Pl (xi,x 2 ) dp i(xi) dp 2 (x 2 ) da 


+ p a- 


p pi i 


' S(f)(xi,X 2 )>a 


S(f) P2 (x i,x 2 ) dpi{xi) dp 2 {x 2 ) da 


C\\f\\ 


Hp{X 1 xX 2 ) 


for all p with p 2 < p < p\. Hence T is bounded from H P (X\ x X 2 ) to L P (X\ x X 2 ), as 
required. 

(b) We turn to the second assertion. For each given A > 0 and / G H p (Xi x X 2 ), by the 
Calderon-Zygmund decomposition again we have 


p({(x!,x 2 ) : jS(Tf)(x 1 ,x 2 )j > a}) 

< p({(x!,x 2 ) : \S(Tg)(x 1 ,x 2 )\ > ^}) + p({(x u x 2 ) : |5(T6)(xi,x 2 )| > |}) 

< Ca- p '\\Tg\\ p ^ 1(XixX2) + Ca 

< Ca~ pi 


P2 \\ Tb \\ P HP 2 ( Xl xX 2 ) 


P hpi( Xi xx 2 ) + Ca P2 \ 


\P2 

\HP2(X 1 xX 2 ) 


< Ca~ Pl 

+ Ca~ P2 


S(f)(xi,x 2 )<a 


2 ) dpi{xi) dp 2 (x 2 ) 


S(f)(xi,x 2 )>a 


S(f) P2 (xi,x 2 ) dpi(xi) dp 2 (x 2 ). 


Therefore \\S(Tf)\\ LP ( Xl xx 2 ) < C\\S(f)\\HP( Xl xx 2 )- Hence \\Tf\\ H p( Xl xx 2 ) < C < \\f\\Hp(x 1 xx 2 ) 
for all p with p 2 < p < pi, as required. □ 
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